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Abstract 

We present a lattice-QCD calculation of the B —)• iriu semileptonic form factors and a new de¬ 
termination of the CKM matrix element \Vub\- We use the MILC asqtad 2+1-flavor lattice configu¬ 
rations at four lattice spacings and light-quark masses down to 1/20 of the physical strange-quark 
mass. We extrapolate the lattice form factors to the continuum using staggered chiral perturbation 
theory in the hard-pion and SU(2) limits. We employ a model-independent z parameterization to 
extrapolate our lattice form factors from large-recoil momentum to the full kinematic range. We 
introduce a new functional method to propagate information from the chiral-continuum extrapola¬ 
tion to the z expansion. We present our results together with a complete systematic error budget, 
including a covariance matrix to enable the combination of our form factors with other lattice-QCD 
and experimental results. To obtain \Vu\,\-, we simultaneously fit the experimental data for the 
B —)• niiy differential decay rate obtained by the BaBar and Belle collaborations together with our 
lattice form-factor results. We find \Vub\ = (3.72 + 0.16) x 10“^ where the error is from the combined 
fit to lattice plus experiments and includes all sources of uncertainty. Our form-factor results bring 
the QCD error on \ Vub\ to the same level as the experimental error. We also provide results for the 
B —)■ niu vector and scalar form factors obtained from the combined lattice and experiment fit, 
which are more precisely-determined than from our lattice-QCD calculation alone. These results 
can be used in other phenomenological applications and to test other approaches to QCD. 

PACS numbers: 13.20.He, 12.38.Gc, 12.15.Hh 


I. INTRODUCTION 

The Cabibbo-Kobayashi-Masakawa (CKM) matrix [1, 2] element \Vub\ is one of the fun¬ 
damental parameters of the Standard Model and is an important input to searches for CP 
violation beyond the Standard Model. Constraints on new physics in the flavor sector are 
commonly cast in terms of over-constraining the apex of the CKM unitarity triangle. In 
contrast to the well-determined angle (3 of the unitarity triangle, the opposite side \Vub/Vcb\ 
is poorly determined, and the uncertainty is currently dominated by \Vub\- This is due to 
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the fact that charmless decays of the B meson have far smaller branching fractions than the 
charmed decays, as well as the fact that the theoretical calculations are less precise than 
for sin2/?, |14s|, or |14b|- Currently the most precise determination of \Vub\ is obtained from 
charmless semileptonic B decays, using exclusive or inclusive methods that rely on the mea¬ 
surements of the branching fractions and the corresponding theoretical inputs. Exclusive 
determinations require knowledge of the form factors, while inclusive determinations rely 
on the operator product expansion, perturbative QCD, and non-perturbative input from 
experiments. There is a long standing discrepancy between \Vub\ determined from inclusive 
and exclusive decays: the central values from these two approaches differ by about 3 ct. It 
was argued in Ref. [3] that this tension is unlikely to be due to new physics effects, and 
it is therefore important to examine the (theoretical and experimental) inputs to the \Vub\ 
determinations. With the result obtained in this paper, the tension is reduced to 2.da. 

In the limit of vanishing lepton mass, the Standard Model prediction for the differential 
decay rate of the exclusive semileptonic B —)■ vrf'z/ decay is given by 


dT{B ^ TTiu) Gl\Vub\ 


dq^ 


247r^ 


|p.ri/+(r)n 


( 1 . 1 ) 


where \Pt^\ = [{^"3 + ~ ~ is the pion momentum in the R-meson 

rest frame. To determine \Vub\-, the form factor |/+(g^)| must be calculated with nonpertur- 
bative methods. The hrst unquenched lattice calculations of |/+(g^)| with 2+1 dynamical 
sea quarks were performed by HPQCD [4] and the Fermilab/MILC collaborations [5] several 
years ago. Here we extend and improve Ref. [5] in several ways. 

The most recent exclusive determination of \ Vub\ from the Heavy Flavor Averaging Group 
(HFAG) [6] is based on combined lattice plus experiment hts and yields | Wb| = (3.28+0.29) x 
10“^, where the error includes both the experimental and theoretical uncertainties. The 
experimental data included in the average are the BaBar untagged six-g^-bin data [7], the 
BaBar untagged twelve-g^-bin data [8], the Belle untagged data [9], and the Belle hadronic 
tagged [10] data. The theoretical errors on the form factors from lattice QGD [5] are currently 
the dominant source of uncertainty in \Vub\ [Hj- Hence a new lattice calculation of 
with improved statistical and systematic errors is desirable To compare, the value of \Vub\ 
from the inclusive method quoted by HFAG is about (4.40 + 0.20) x 10“^ [6] using the theory 
of Ref. [15]. 


Note that there are several other efforts with 2 [12] and 2+1 flavors of sea quarks [13, 14[. 
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In this paper, we present a new lattice-QCD calculation of the B —)■ 'K^l> semileptonic 
form factors and a determination of \Vub\- Our calculation shares some features with the 
previous Fermilab/MILC calculation [5] but makes several improvements. We quadruple 
the statistics on the previously used ensembles and improve our strategy for extracting the 
form factors by including excited states in our three-point correlator analysis. In addition, 
we include twice as many ensembles in this analysis. The new ensembles have smaller 
lattice spacings, with the smallest lattice spacing decreased by half. This analysis also 
includes ensembles with light sea-quark masses that are much closer to their physical values 
{mi/ms = 0.05 versus 0.1). The smaller lattice spacings and light-quark masses provide much 
better control over the dominant systematic error due to the chiral-continuum extrapolation. 
We find that heavy-meson rooted staggered chiral perturbation theory (HMrSyPT) in the 
SU(2) and hard-pion limits provides a satisfactory description of our data. All together, 
these improvements reduce the error on the form factors by a factor of about 3. Finally, we 
introduce a new functional method for the extrapolation over the full kinematic range. 

The determination of \Vub\ from a combined fit to our lattice form factors together with 
experimental measurements also yields a very precise determination of the vector and scalar 
form factors over the entire kinematic range. These form factors will be valuable input to 
other phenomenological applications in the Standard Model and beyond. An example is the 
rare decay B —)• , which we will discuss in a separate paper. 

Because our primary goal in this work was a reliable and precise determination of \Vuh\, 
we employed a blinding procedure to minimize subjective bias. At the stage of matching 
between the lattice and continuum vector currents, a slight multiplicative offset was applied 
to the data that was only known to two of the authors. The numerical value of the blinding 
factor was only disclosed after the analysis and error-estimation procedure, including the 
determination of \Vub\, were essentially finalized. 

This paper is organized as follows. In Sec. II, we present our calculation of the form fac¬ 
tors. We describe the lattice actions, currents, simulation parameters, correlation functions 
and fits to extract the matrix elements, renormalization of the currents, and adjustment of 
the form factors to correct for quark-mass mistunings. In Sec. Ill, we present the combined 
chiral-continuum extrapolation, followed by an itemized presentation of our complete error 
budget in Sec. IV. We then extrapolate the form factors to the full range through the 
functional z expansion method in Sec. V. We also perform fits to lattice and experimental 


5 


data simultaneously, to obtain \Vub\- We conclude with a comparison to other results and 
discussion of the future outlook in Sec. VI. Preliminary reports of this work can be found in 
Refs. [16, 17]. 

II. LATTICE-QCD SIMULATION 

In this section, we describe the details of the lattice simulation. We briefly describe the 
calculation of the form factors in Sec. II A. We also calculate the tensor form factor, which 
follows a analysis similar to that of the vector and scalar form factors. The tensor form 
factor enters the Standard-Model rate for B —)■ decay, and our hnal result for Jt will 

be presented in a forthcoming paper. In Sec. IIB, we introduce the actions and simulation 
parameters used in this analysis. This is followed, in Sec. IIC, by a brief discussion of the 
currents and lattice correlation functions. The correlator hts to extract the lattice form 
factors are provided in Sec. IID. In Sec. HE, we discuss the renormalization of the lattice 
currents. In Sec. IIF, we correct the form factors a posteriori to account for the mistuning 
of the simulated heavy 6-quark mass. 

A. Form-factor definitions 

The vector and tensor hadronic matrix elements relevant for I? —)■ tt semileptonic decays 
can be parameterized by the following three form factors: 



( 2 . 1 ) 


( 2 . 2 ) 


where and = iqa^'^h. In lattice gauge theory and in chiral perturbation 

theory, it is convenient to parameterize the vector-current matrix elements by [18] 


{ir(p,)|V'‘|B(PB)) = [v‘‘hiE,)+p:jAE,)] , 


(2.3) 


where = p^/Ms is the four velocity of the B meson and p^_|_ = p^ — ■ v)v^ is the 

projection of the pion momentum in the direction perpendicular to v^. The pion energy is 
related to the lepton momentum transfer by • u = (M^ -|- — (f)/ (2Mb). With 
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this setup, we have 


fwiE^) 

f±{En) 

Mq^) 


{7r{p^)\V^\B{ps)) 
{n{p^)\V^\B{pB)) ^ 

pV 

Mb + M^ {7r{p^)\r^^\B{pB)) 1 
a/2Mb \/2Mb P\: 


(2.4) 

(2.5) 

( 2 . 6 ) 


where no summation is implied by the repeated indices here. The form factors /+ and /o 


are 


/+(«") = [ft + (Mb - E,)B] , ( 2 . 7 ) 

+ (El - MDB] . ( 2 . 8 ) 


B. Actions and parameters 

The lattice gauge-field conhgurations we use have been generated by the MILC Collab¬ 
oration [19-21], and some of their properties are listed in Table I. These twelve ensembles 
have four different lattice spacings ranging from a ~ 0.12 fm to a ~ 0.045 fm with several 
light sea-quark masses at most lattice spacings in the range 0.05 < am'Jam'f^ < 0.4. The 
parameter range is shown in Fig. 1. We use the Symanzik-improved gauge action [22-24] for 
the gluons and the tadpole-improved (asqtad) staggered action [25-30] for the 2-1-1 flavors of 
dynamical sea quarks and for the light valence quarks. Both Table I and Fig. 1 also indicate 
the ensembles used in the previous Fermilab/MILC calculation [5]. The current analysis 
benehts from an almost quadrupled increase in the statistics over that of Ref. [5] , as well as 
hner lattice spacings and lighter sea-quark masses. All ensembles have large enough spatial 
volume, MjrL > 3.8, such that the systematic error due to hnite-size effects is negligible 
compared to other uncertainties. 

In this calculation, we work in the full-QCD limit, so that the light valence-quark masses 
ami are the same as the light sea-quark masses am[, which are degenerate. For the bottom 
quarks, we use the Fermilab interpretation [31] of the Sheikholeslami-Wohlert clover action 
[32]. In Table II, we list parameters for the valence quarks. 

Table III lists the values of ri/a on each ensemble, along with other derived parameters, 
where ri is the characteristic distance between two static quarks such that the force between 
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Table I. Parameters of the MILC asqtad gauge-field ensembles used in this analysis. From left 
to right: approximate lattice spacing a in fm, the (light/strange)-quark mass ratio am'Jam'y^, the 
coupling constant /3, the tadpole parameter uq determined from the plaquette, lattice volume, the 
number of configurations iVcfg, M-j^L (L is the spatial length of the lattice), and the number of 
configurations of the four ensembles that were used in Ref. [5]. 


Ria(fm) 

am'i/am'^ 

13 

Uo 

volume 

Ncfg 

M^L 

iVrfg(Ref. [5])) 

0.12 

0.01/0.05 

6.760 

0.8677 

20^ X 64 

2259 

4.5 

592 


0.007/0.05 

6.760 

0.8678 

20^ X 64 

2110 

3.8 

836 


0.005/0.05 

6.760 

0.8678 

243 X 64 

2099 

3.8 

529 

0.09 

0.0062/0.031 

7.090 

0.8782 

28^ X 96 

1931 

4.1 

557 


0.00465/0.031 7.085 

0.8781 

32^ X 96 

984 

4.1 



0.0031/0.031 

7.080 

0.8779 

40^ X 96 

1015 

4.2 



0.00155/0.031 7.075 0.877805 

64^ X 96 

791 

4.8 


0.06 

0.0072/0.018 

7.480 

0.8881 

48^ X 144 

593 

6.3 



0.0036/0.018 

7.470 

0.88788 

48^ X 144 

673 

4.5 



0.0025/0.018 

7.465 

0.88776 

56^ X 144 

801 

4.4 



0.0018/0.018 

7.460 

0.88764 

64^ X 144 

827 

4.3 


0.045 

0.0028/0.014 

7.810 

0.89511 

64^ X 192 

801 

4.6 



them satisfies r^F(ri) = 1.0 [33, 34]. The absolute lattice scale ri is obtained by comparing 
the Particle Data Group value of /jr with lattice calculations of ri/^r from MILC [35] and 
HPQCD [36], yielding the absolute scale ri = 0.3117(22) fm [37]. This value is consistent 
with the independent, but less precise, determination ri = 0.323(9) from RBC/UKQCD 
using domain-wall fermions [38]. 

C. Currents and correlation functions 

We calculate the two-point and three-point functions 

Gp(f;p) = ^e*P-*(Op(0,0)Oj,(f,a;)), and (2.9) 

X 

Cj{t,T;p) = ^e-»'“{O,(0,0)J(t.y)Ot(r,ii:)>. 

x,y 

8 


( 2 . 10 ) 
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Figure 1. Lattice spacings and light-quark masses used in this analysis. The area of each hlled 
disk is proportional to the number of conhgurations in the ensemble. Open black circles indicate 
ensembles use in the analysis of Ref. [5]; those with (a, m;/m^) = (0.12 fm, 0.4) and (0.09 fm, 0.4) 
used in Ref. [5] (open circles without disks) are not used in this analysis. 

where P = B, labels the pseudoscalar meson, the operators Op (Op) annihilate (create) 
the states with the quantum numbers of the pseudoscalar meson P on the lattice, and 
J = are the lattice currents. 

For the B meson, we use a mixed-action interpolating operator Op which is a combination 
of a Wilson clover bottom quark and a staggered light quark [5]; 

C>b(x) = J2y^iy)S{y,xh5^{x)x{x), (2.11) 

where L2(a;) = 71 ^ 72 ^ 73 ^ 74 "^, x = and S{x,y) is a smearing function. For the pion, we 

use the operator 

= {-^)^*=^'''x{x)x{x), ( 2 . 12 ) 

which is constructed from two 1 -component staggered quarks. 

The current operators are constructed in a similar way: 

V^{x) = T(x) 7 '^r 2 (a:)x(a;), and 
= T(x)a^'^f2(a;)x(a;), 
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(2.13) 

(2.14) 



Table II. Heavy-quark masses and other parameters used in the simulation. Starting in the third 
column: the clover parameter cswj the simulation 6-quark mass parameter the current rotation 
parameter the number of sources Ngj-c and the two source-sink separations T. Note that we use 
the same valence light-quark mass as in the sea except the a = 0.09fm, = 0.00465/0.031 

ensemble where a slightly different valence mass ami = 0.0047 is used. 


« a(fm) 

am'Jam'^ 

csw 

d[ 

Asrc 

T 

0.12 

0.01/0.05 

1.531 

0.0901 0.09334 

4 

18,19 


0.007/0.05 

1.530 

0.0901 0.09332 

4 

18,19 


0.005/0.05 

1.530 

0.0901 0.09332 

4 

18,19 

0.09 

0.0062/0.031 

1.476 

0.0979 0.09677 

4 

25,26 


0.00465/0.031 

1.477 

0.0977 0.09671 

4 

25,26 


0.0031/0.031 

1.478 

0.0976 0.09669 

4 

25,26 


0.00155/0.031 

1.4784 

0.0976 0.09669 

4 

25,26 

0.06 

0.0072/0.018 

1.4276 

0.1048 0.09636 

4 

36,37 


0.0036/0.018 

1.4287 

0.1052 0.09631 

4 

36,37 


0.0025/0.018 

1.4293 

0.1052 0.09633 

4 

36,37 


0.0018/0.018 

1.4298 

0.1052 0.09635 

4 

36,37 

0.045 

0.0028/0.014 

1.3943 

0.1143 0.08864 

4 

48,49 


where the heavy quark field spinor H/ is rotated to remove tree-level 0{a) discretization 
effects, via [31] 


d'(a;) = (1 -k acii7 • Aat)V’- (2-15) 

Figure 2 illustrates the three-point correlation function used to obtain the lattice form 
factors. The current operator J is inserted between the b- and /-quark lines. The three- 
point functions depend on both the current insertion time t and the temporal separation 
T between the tt and B mesons. The signal to noise ratio is largely determined by T. A 
convenient approach is to fix the source-sink separation T in the simulations and then insert 
the current operators at every time slice in between. The source-sink separations T at differ¬ 
ent lattice spacings, sea-quark masses, and recoil momenta are chosen to be approximately 
the same in physical units. To minimize statistical uncertainties and reduce excited-state 
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Table III. Derived parameters from the simulation. Starting in the third column: relative scale 
ri/a, the Goldstone pion mass root-mean-square (RMS) pion mass and the critical 

hopping parameter Kent which enters our definition of the heavy-quark mass. 


« o(fm) 


rtja 

M^(MeV) M^^S(MeV) 

^crit 

0.12 

0.01/0.05 

2.7386 

389 

532 

0.14091 


0.007/0.05 

2.7386 

327 

488 

0.14095 


0.005/0.05 

2.7386 

277 

456 

0.14096 

0.09 

0.0062/0.031 

3.7887 

354 

413 

0.139119 


0.00465/0.031 3.7716 

307 

374 

0.139134 


0.0031/0.031 

3.7546 

249 

329 

0.139173 


0.00155/0.031 3.7376 

177 

277 

0.139190 

0.06 

0.0072/0.018 

5.3991 

450 

466 

0.137582 


0.0036/0.018 

5.3531 

316 

340 

0.137632 


0.0025/0.018 

5.3302 

264 

291 

0.137667 


0.0018/0.018 

5.3073 

224 

255 

0.137678 

0.045 

0.0028/0.014 

7.2082 

324 

331 

0.136640 



Figure 2. Illustrative diagram for the three-point correlation functions. 

contamination, we tested data with different source-sink separations before choosing those 
shown in Table II. The B meson is at rest in our simulation, while the daughter pion is 
either at rest or has a small three-momentum. The light-quark propagator is computed 
from a point source so that one inversion of the Dirac operator can be used to obtain mul- 
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tiple momenta. The spatial source location is varied randomly from one configuration to 
the next to minimize autocorrelations. The 6-quark source is always implemented with 
smearing based on a Richardson IS wave function [39] after fixing to Coulomb gauge. We 
compute both the two-point function C.^{t]p) and three-point function Cj{t,T]p) at several 
of the lowest possible pion momenta in a finite box: p = (27r/L)(0, 0, 0), (27r/L)(l, 0, 0), 
(27r/L)(l, 1, 0), (27r/L)(l, 1,1), and (27r/L)(2, 0, 0), where contributions from each momen¬ 
tum are averaged over permutations of components. We find the correlation functions with 
momentum (27r/L)(2,0, 0) too noisy to be useful, so we exclude these data from our analysis. 


D. Two-point and three-point correlator fits 


In this subsection, we describe how to extract the desired matrix element from two- and 
three-point correlation functions. With our choice for the valence-quark actions and for the 
interpolating operators, the two- and three-point functions take the form [40] 




n=0 

oo 




(2.16) 


Cj{t,T;p) = T-0^ 

m,n=0 

(2.17) 


where W is the temporal length of the lattice and 

|(0|Op|PW(p))| 


-P (P) = 


M 


(mn) 

J 


2Ejr\p) 

2Ei^\p)^/2M^^ 


(2.18) 

(2.19) 


Note that due to the staggered action used for the light quarks, the meson interpolating 
operators also couple to the positive parity (scalar) states which oscillate in Euclidean times 
t and T with the factors (—l)”h+i) and 

Our goal is to extract the ground state matrix element from these correlation 

functions. To suppress the contributions from the positive parity states to the ratio, we follow 
the averaging procedure of Ref. [5] , which exploits the oscillating sign in their Euclidean time 
dependence. The time averages can be thought of as a smearing over neighboring time slices 
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{t, t + 1, t + 2} X {T, T + 1} to significantly reduce the overlap with opposite-parity states. 
Denoting the averaged correlators by Cp and Cj, we then use the ratio [5] 



( 2 . 20 ) 


where E^\p) and are the ground-state pion energy and i?-nieson rest mass, respec¬ 


tively. The uncertainty in the i?-meson rest mass has signihcant impact on the ratio i?j, 
so we follow a two-step procedure. We hrst determine the pion and i?-meson ground-state 
energy as precisely as possible using the corresponding two-point functions. We then feed 
these ground-state energies into the ratio Rj, preserving the correlations with jackknife 
resampling. 

For the pion two-point functions at zero momentum, the oscillating states — the terms 
in Eq. (2.16) with odd powers of (—1) — do not appear. Thus, we fit the pion two-point 
functions using Eq. (2.16) with the lowest two non-oscillating states (n = 0, 2). For the two- 
point functions with nonzero momentum, the contribution from oscillating states is small 
but noticeable. We hud that we only need to include the lowest three states {n = 0,1, 2) 
in the hts. Because the momenta we consider are typically small compared to 27r/a, the 
continuum dispersion relation is satished within statistical errors, as shown in Fig. 3. In the 
main analysis, we therefore use the mass from the zero-momentum ht and the continuum 
dispersion relation to set (p) = -I- for non-zero momentum. Because the zero- 

momentum energy has signihcantly smaller statistical error than that of nonzero momentum, 
using this choice and the dispersion relation for nonzero-momentum energy leads to a more 
stable and precise determination of Table IV lists the relevant fit ranges for the 

two-point hts. In the two-point correlators (except the zero-momentum pion two-point 
correlators), the noise grows rapidly with increasing t, the distance away from the pion 
source in the temporal direction. The data points at large t are not useful, and including 
them would lead to a larger covariance matrix which would be difficult to resolve given the 
limited number of conhgurations. We choose the upper end of the ht ranges fmax such that 
the relative error does not exceed 20%. The lower end fmin is chosen such that the excited 
state contamination is sufficiently small, he., the resulting central values of the ground state 
energy are stable against variations in tmin as shown in Fig. 4 (left). 

In our analysis, there are two places where quantities from the S-meson two-point func- 
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C\J 

LU 



3 

O 

o 

3 

O 

o 



(|p|L/(27t))^ 


(|p|L/(27t)r 


Figure 3. Dispersion-relation tests on the ensemble a ~ 0.12 fm, = 0.1. Left: fit results of 

from two-point functions are compared with the continuum dispersion relation |pp -|- M^. 
Right: wave-function overlaps (0|O.n-|7r(p)) at different momenta are compared. The dotted lines 
show the expected size of the leading momentum-dependent discretization errors based on power 
counting, which are of 0{a‘^\p\‘^a‘^). 


tions are needed. The first is for in Eq. (2.20). The second is for the 5-nieson excited 
state energy in Eq. (2.21) below. For the determination of in Eq. (2.20) we use two- 
point functions constructed with a IS-smearing function in the interpolating operators for 
the source and sink. The IS-smeared operator has good overlap with the ground state and 
a much smaller overlap with the excited states than the local source operator, thus reducing 
excited-state contributions to the corresponding correlators. We fit the IS-IS smeared B- 
meson two-point correlators with relatively large tmin to only two states (n < 2 in Eq. (2.16)). 
To choose fmax, we again apply the 20%-rule on the relative error. The lower bound fmin 
is chosen in a manner similar to the pion two-point fits and the stability plot is shown in 
Fig. 4 (right). The chosen fit ranges are shown in Table IV. 

We test for autocorrelations by blocking the configurations on each ensemble with dif¬ 
ferent block sizes, and then using a single-elimination jackknife procedure to propagate the 
statistical error to the two-point correlator fits for and M^\ We do not observe any no¬ 
ticeable autocorrelations in all the ensembles we use, as illustrated in Fig. 5 for the coarsest 
and finest ensembles, and choose not to block the data. 

The ratios in Eq. (2.20) have the advantage that the wavefunction overlap factors Zp 
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Figure 4. Pion (left) and B meson (right) mass in lattice units from the two-point correlator fits 
versus different choices of tmin for the a ~ 0.12 fm, = 0.1 ensemble. The tmin marked in red 

(filled) is that of our preferred fit. The pion and B meson correlators are fit to 2+0 and 1+1 states 
with tmax = 30, 21, respectively. 


Table IV. Fit ranges [tminjimax] of the two-point correlator fits used to obtain the rest masses of 
the pion and B mesons. 


« a(fm) 

1 am'Jam'j^ 

C-K (2+0) 

(1+1) 

0.12 

0.01/0.05 

[6, 30] 

[9, 21] 


0.007/0.05 

[6, 30] 

[9, 21] 


0.005/0.05 

[6, 30] 

[9, 21] 

0.09 

0.0062/0.031 

[9, 47] 

[12, 32] 


0.00465/0.031 

[9, 47] 

[12, 29] 


0.0031/0.031 

[9, 47] 

[13, 29] 


0.00155/0.031 

[9, 47] 

[14, 29] 

0.06 

0.0072/0.018 

[13, 71] 

[14, 41] 


0.0036/0.018 

[13, 71] 

[14, 41] 


0.0025/0.018 

[13, 71] 

[14, 41] 


0.0018/0.018 

[13, 71] 

[15, 41] 

0.045 

0.0028/0.014 

[17, 74] 

[17, 61] 


15 





























0.1602 


2 


0.1601 

0.16 

0.1599 

0.1598 

0.1597 

0.1596 


fitted I— e 

error{blocking)/error(non-blocking) . 


c 

c 

) . 

> 

= . 

C) 







0 1 2 3 4 5 6 7 

Block size 


2 


1 


0 


0.0714 
0.0712 

0.071 

0.0708 

0.0706 

0 1 2 3 4 5 6 7 

Block size 


fitted I— 

error{blocking)/error(non-blocking) 





O .-.-.C 

T .' 

o Q c 

) « 



1 


1.845 

1.84 

1.835 

1.83 

1.825 

0 1 2 3 4 5 6 7 

Block size 


' fitted Mg I— '-Q- 


error(blocking)/error{non-blocking) 


cj . o . o 





« . c 


) “ 


1 


0.934 
0.9335 
0.933 
0.9325 
0.932 
0.9315 
0.931 

0 1 2 3 4 5 6 7 

Block size 


fitted Mg I — e- 


error{blocking)/error(non-blocking) 


o —xr 


2 


1 


0 


Figure 5. Fitted and Mb in lattice units versus the block size for the a ~ 0.12 fm, = 0.1 

ensemble (top) and for the a ~ 0.045 fm, = 0.2 ensemble (bottom). The dashed line shows 

the ratio of the ht errors from the blocking and non-blocking data, which can be read off from the 
right y-axis. 


cancel, but the trade-off is that we need an additional factor — the square root term on the 
right-hand side — to remove the leading t dependence in the ratio. If the lowest lying states 
dominated the ratio Rj, then it would be constant in t and proportional to the lattice form 
factor fj. The subscript J now runs over T, ||, and T, corresponding to the operators Id*, 
V‘^, and T"^*, respectively. Our previous analysis employed a simple plateau fit constant in 
time. With our improved statistics, the small excited-state contributions to the ratio are 
significant and cannot be neglected. On the other hand, even with our improved statistics, 
we find that contributions to Rj from wrong-parity states are still negligible. We use two 
different fit strategies to remove excited state contributions and use the consistency between 
them as an added check that any remaining excited state contamination is negligibly small. 

The first strategy starts with the ratio in Eq. (2.20) and minimally extends the plateau 
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Figure 6. Ground-state energy (bottom) and the first excited-state splitting (top) of the B meson 
on the a ~ 0.12 fm, = 0.1 ensemble, from different fits to local, IS-smeared and both 

two-point correlators. 

fitting scheme by including the first excited state of the B meson in the following form; 

= ff h + , (2.21) 

hj 

where Aj and /j"* are unconstrained £t parameters, AMb = is the lowest 

energy splitting of the pseudoscalar B meson, and the prefactors are h\\ = 1, h± = and 
hx = {y/2MBp\^) / {Mb + MA)- We choose the £t ranges for Rj such that contributions from 
pion excited states to Rj can be neglected. The ht parameter AMb is determined by the 
B-meson two-point correlators. In practice, we fit the ratio in Eq. (2.21) along with the 
i?-meson two-point correlation functions with AMb as a common parameter. We find it 
benehcial in the combined £t to include both the local and smeared two-point correlation 
functions. We use 2-1-2 states for both correlators, but use a different set of fit ranges (listed 
in Table V). The results of these two-point hts are shown in Fig. 6. The agreement in the 
i?-meson energies between the separate and combined hts is very good, but the combined 
ht leads to smaller errors. 

To summarize our strategy, for the case of zero momentum, we ht the ratio i?|| (f) together 
with the local and smeared 5-meson two-point correlators simultaneously. For 

non-zero momentum p, AMb is common to all three ratios, R\\, R±, Rt- Thus, we perform 
a combined ht to the hve quantities: i?||, R±, Rt, and Figure 7 shows an 
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Table V. The fit ranges [tminj^max] of the combined two-point and three-point ratio fits to obtain 
the lattice form factors. 


« a(fm) 

1 am'Jam!^ 

'-'B 

'-'B 

5|| 

R± 

Rt 

0.12 

0.01/0.05 

[9, 23] 

[7, 21] 

[6,12] 

[6,12] 

[6,12] 


0.007/0.05 

[9, 23] 

[7, 21] 

[6,12] 

[6,12] 

[6,12] 


0.005/0.05 

[9, 23] 

[7, 21] 

[6,12] 

[6,12] 

[6,12] 

0.09 

0.0062/0.031 

[12, 32] 

[9, 32] 

[9,16] 

[9,16] 

[9,16] 


0.00465/0.031 

[12, 32] 

[9, 29] 

[9,16] 

[9,16] 

[9,16] 


0.0031/0.031 

[12, 32] 

[9, 29] 

[9,16] 

[9,16] 

[9,16] 


0.00155/0.031 

[12, 29] 

[9, 29] 

[9,15] 

[9,15] 

[9,15] 

0.06 

0.0072/0.018 

[13, 41] 

[9, 41] 

[12,22] 

[12,22] 

[12,22] 


0.0036/0.018 

[13, 41] 

[9, 41] 

[12,22] 

[12,22] 

[12,22] 


0.0025/0.018 

[13, 41] 

[9, 41] 

[12,22] 

[12,22] 

[12,22] 


0.0018/0.018 

[13, 41] 

[9, 41] 

[12,21] 

[12,21] 

[12,21] 

0.045 

0.0028/0.014 

[16, 61] 

[10, 61] 

[16,26] 

[16,26] 

[16,26] 


example of these fits. Figure 8 shows the stability plots of against the variations in the 
fit ranges of the ratio fits, and the variations in the fit ranges of both two-point correlators. 
The preferred fit ranges are set to be in the stable region upon these variations. 

Our second fit strategy includes excited-state contributions from both the pion and the 
B meson. It starts with a different ratio, without time averages, which ensures that there 
are enough data points to constrain all the parameters: 


Rj{t) 


Cj{t,T-p) 

c^{t-p)CB{T -t-oy 


( 2 . 22 ) 


where the two- and three-point correlators are defined in Eqs. (2.9) and (2.10). We fit Rj{t) 
with all the possible states with m,n < 2 in Eqs. (2.16) and (2.17), combining the fits to 
the pion and 5-meson two-point correlators. We compare the fit results of the two different 
fit schemes in Fig. 9. The first (simple) fit model described in Eq. (2.21) gives, fitting either 
simultaneously or individually to the three lattice form factors fj, results that are consistent 
with the second fit model that includes the full set of first excited states in Eq. (2.22). In 
contrast, the plateau fits to Rj defined in Eq. (2.20) yield results that are as much as one 


18 












5 p=^(0,0,0), p value=0.82 
5 p=|L(i,o,0), p value=0.29 (combined fit) 
(J p=^(l,l,,0), p vaiue=0.98 (combined fit) 
^ p=^(l,l,l), p vaiue=0.99 (combined fit) 




Figure 7. Fits of ratios R±/h± (top left), R\\/h^^ (bottom left) and Rt/Ht (bottom right) for 
the a ~ 0.12 fm, = 0.1 ensemble. The open data points with error bars are the ratios 

constructed from two- and three-point correlators with various momenta. The fit results are shown 
as solid lines with dashed error widths. The shaded bands indicate the range of data employed in 
the fit. The resulting lattice form factor fj are shown by the colored bars on the left of the plots. 

statistical a smaller. In summary, we find that the first fit strategy described by Eq. (2.21) 
is sufficient to remove contributions from excited states, and we therefore adopt this method 
for the main analysis. 


E. Matching 

We match the lattice currents to continuum QCD with the relation, 

J = Zj„J, (2.23) 
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Figure 8. Left: fitted ratio R±/h± on the a ~ 0.12 fm, m^/m^ = 0.1 ensemble with momentum 
p = (27r/L)(l, 0, 0) for different combinations of tminjimax of the three-point correlator in the 
combined fit, where circles are results of varying fmin with fixed fmax = 12 and squares are that 
of varying fmax with fixed fmin = 6. Right: the same R±/h± upon variations of fmin of the two 
i3-meson two-point correlators (local and IS-smeared) in the combined fit, where tmax is given in 
Table V. The filled points in both plots show the preferred fit ranges and the dashed lines indicate 
the p values of the fits. 



Figure 9. Fitted lattice form factors fj at different momenta: from top to bottom, (27r/L)(0, 0, 0) 
(/|| only), (27r/L)(l, 0, 0), (27r/L)(l, 1,0) and (27r/L)(l, 1,1). Each compares the following differ¬ 
ent fit schemes: separate (squares) or combined (circles) fits of the three lattice form factors fj 
using Eq. (2.21), fit including the full first excited states as defined in Eq. (2.22) (upward-pointing 
triangles) and the simple plateau fit (downward-pointing triangles). 
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where J and J denote the vector or tensor currents in the continuum and lattice theories, 
respectively, and “=” means “has the same matrix elements” [41], We calculate the current 
renormalization with the mostly nonperturbative renormalization method [18, 42], 

(2-24) 

where ZyA and ZyA are the matching factors for the corresponding flavor-conserving vector 
currents. These factors capture most of the current renormalization. The remaining flavor 
off-diagonal contribution to the matching factor, is close to unity. 

We calculate the factors ZyA^ and ZyA nonperturbatively for each ensemble by computing 
the matrix elements of the flavor-conserving vector currents and using the relations 

l = Zypir\V‘\7r), (2.26) 

l = ZypB.\V*\B.), (2,26) 

where the lattice current is a bilinear of light staggered quark helds and is a bilinear 
of clover heavy quark helds. The factors ZyA and ZyA are listed in Table VI. Because there 
is very little m; dependence in the factor ZyA, we use the same ZyA for ensembles with 
different light quark masses but the same lattice spacing. The factor ZyA^ depends crucially 
on the heavy b quark mass, though it has negligible light quark mass dependence. 

We use lattice perturbation theory [43] to compute the remaining renormalization factors 
Pj at one-loop. Due to the cancellation of the tadpole contributions in the radiative correc¬ 
tions to the left and right side of Eq. (2.24), the factors pj are very close to one. They have 
the perturbative expansion 

pj = l + ay{q*)p^]^ + 0{al), (2.27) 

where we take the strong coupling in the V-scheme [43] at a scale q* that corresponds to 
the typical gluon loop momentum. In practice, we choose q* = 2/a. The details of the 
calculation of the one-loop coefficients will be presented elsewhere. The values used in 
this work are shown in Table VI. 

F. Heavy-quark mass correction 

In the clover action, the hopping parameter corresponds to the bare 6-quark mass. 
When we started generating data for this analysis, we had a good estimate for the bottom- 
quark on each ensemble, but not the final tuned values, which were obtained as described 
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Table VI. The parameters for the renormalization of the form factors. The first two columns label 
the ensemble with its approximate lattice spacing and the sea light- and strange-quark mass ratio. 
The third column is the simulation The fourth and hfth columns are the nonperturbative heavy- 
heavy and light-light renormalization factors Zyi^, Zy4. The sixth, seventh, and eighth columns are 
the one-loop estimates of py4, pyi and pT, respectively. The tensor current has a nonzero anomalous 
dimension; the numbers reported here match to the MS scheme at renormalization scale p = m 2 . 


which corresponds to the pole mass. Note that with our convention Zy4 and Zy4^ are normalized 
so that at tree-level Zyl = 2 and Zy\ = 1 — 6ttoK^ [37]. As a result, Zy bb only approaches 1 in the 

* LL ^bb ^ 

limit rrih —)• 00 . 


a(fm) 

am'Jam'^ 

4 



PvA 

Pvii 

PTbi 


0.12 0.01/0.05 0.0901 0.5015(8) 1.741(3) 1.006214 0.973023 1.033350 

0.007/0.05 0.0901 0.5015(8) 1.741(3) 1.006252 0.973109 1.033280 

0.005/0.05 0.0901 0.5015(8) 1.741(3) 1.006197 0.973082 1.033270 

0.01/0.05 0.0860 0.5246(9) 1.741(3) 1.012999 0.977290 1.030650 

0.01/0.05 0.0820 0.5469(10) 1.741(3) 1.018261 0.980129 1.028960 


0.09 0.0062/0.031 

0.00465/0.031 
0.0031/0.031 
0.00155/0.031 


0.0979 0.4519(15) 1.776(5) 0.999308 0.975822 1.036590 
0.0977 0.4530(15) 1.776(5) 0.999405 0.975775 1.036390 
0.0976 0.4536(15) 1.776(5) 0.999441 0.975744 1.036350 
0.0976 0.4536(15) 1.776(5) 0.999416 0.975703 1.036390 


0.1048 0.4089(21) 1.807(7) 0.995605 0.979279 1.042390 
0.1052 0.4065(21) 1.807(7) 0.995371 0.979260 1.043160 
0.1052 0.4065(21) 1.807(7) 0.995350 0.979217 1.043190 
0.1052 0.4065(21) 1.807(7) 0.995327 0.979176 1.043250 
0.045 0.0028/0.014 0.1143 0.3564(65) 1.841(6) 0.994195 0.984351 1.058790 


0.06 0.0072/0.018 

0.0036/0.018 
0.0025/0.018 
0.0018/0.018 


in Appendix C of Ref. [44]. We therefore need to adjust the form factors a posteriori to 
account for the slightly mistuned values of Kb- 

The Kb parameters are adjusted so that the corresponding Bg kinetic masses match the 
experimentally-measured value [44] . Table VII shows both the simulation and hnal tuned Kb 
values. For some ensembles, the difference between the two is as large as 7(T of the statistical 
uncertainty associated with the tuning procedure. We study the Kb-dependence of the lattice 
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Table VII. Parameters needed to apply heavy-quark mass corrections. The third column contains 
the value used for the calculation, the fourth column contains the tuned value Ki, with its 
statistical error. Subsequent columns contain the percentage shift in m2 and each of the form 
factors. 


a(Rifm) 

am(/am)j 

K'b Kb 


^(%) 

f(%) 

^(%) 

0.12 

0.01/0.05 

0.0901 0.0868(9) 

10.9 

1.79 

1.55 

1.60 


0.007/0.05 

0.0901 0.0868(9) 

10.9 

1.80 

1.57 

1.58 


0.005/0.05 

0.0901 0.0868(9) 

10.9 

1.81 

1.58 

1.56 

0.09 

0.0062/0.031 

0.0979 0.0967(8) 

4.3 

0.69 

0.60 

0.62 


0.00465/0.031 0.0977 0.0966(8) 

3.9 

0.63 

0.55 

0.56 


0.0031/0.031 

0.0976 0.0965(8) 

3.9 

0.64 

0.56 

0.55 


0.00155/0.031 0.0976 0.0964(8) 

4.2 

0.70 

0.61 

0.59 

0.06 

0.0072/0.018 

0.1048 0.1054(5) 

-2.4 

-0.37 

-0.36 

-0.44 


0.0036/0.018 

0.1052 0.1052(5) 

0.0 

0.0 

0.0 

0.0 


0.0025/0.018 

0.1052 0.1051(5) 

0.4 

0.06 

0.06 

0.06 


0.0018/0.018 

0.1052 0.1050(5) 

0.8 

0.13 

0.11 

0.11 

0.045 

0.0028/0.014 

0.1143 0.1116(4) 

14.3 

2.34 

2.03 

2.10 


form factors by generating data on the a ~ 0.12 fm, = 0.2 ensemble, with two 

additional valnes, 0.0860 and 0.0820, and all other simnlation parameters unchanged. To 
generalize the dependence from this ensemble to others, we work with the quark kinetic 
mass m 2 instead of Kb itself. We expand the form factor f {f = /±,||,r) in about a 
reference point (which corresponds to the tuned Kb) as follows 


/(m2 \ mi, E^) k. /(m^\ mi, E^) -k 
= f{ffi2^,a?,mi,E^) 

= f{fh2^,a^,mi,E^) 1 - 


dfjm^ ,a^,mi,E^) 
dfru^ 


m, 


-1 


m 


-1 


) 


1 + (m-^ 


f dm 2 
dlnf m2 


mn 




d In m 2 m 2 




(2.28) 


where the masses and E^^ are all in ri units. To obtain / at the reference point, we need to 
find the dimensionless normalized slope — (9 In//c? In m2). 

We use exactly the same procedure as described in Sec. IID for k'^, = 0.0901 to obtain the 


23 






















Figure 10. The normalized slopes —(91n//91n?fi2) of the three different form factors /||,_l,t 
several momenta for i? —)■ vr (left) and B ^ K (right) semileptonic decays. The horizontal shaded 
error bands in each plot are the results of correlated hts to all momenta for each form factor. 


B —>■ Tiiu form factors /||,±,t for fhe additional values = 0.0860 and 0.0820. We apply the 
matching factors given in Table VI. Finally, we take m2 to be the kinetic mass corresponding 
to Kb = 0.0868 (the tuned kappa given in Table VII) and use it as the reference point. We 
ht each form factor at each momentum for the three data points to the linear form given in 
Eq. (2.28), taking /(m^^) and —(91n//91nm2) as £t parameters. The result is shown in 
Fig. 10 (left). As shown in the plot, the normalized slope — {din f /din m 2 ) has a very mild 
dependence. Therefore, for each form factor we perform a correlated fit to all momenta 
to obtain a single common normalized slope. The result is shown in Table VIII. Fitting the 
data to a linear form in results in a slope statistically consistent with zero. 

To examine the light-quark mass dependence of the normalized slopes, we repeat the same 
procedure for the B ^ K semileptonic form factors with a heavier daughter valence quark 
amg = 0.0349, which is close to the physical strange-quark mass. The results are plotted in 
Fig. 10 (right). We £t the points of each form factor to a constant and tabulate the results 
in Table VIII. Comparing the normalized slopes for and taking into account 

statistical correlations, we observe a mild but statistically-signihcant light daughter-quark 
mass dependence. So we £t the slopes for and simultaneously to a linear form. 


dlnf mi 

——— — c + a —, 
a In m 2 mg 


(2.29) 


where mi/mg = 0.2 and 1.0 for and respectively. The results for the parameters 

c and d are given in Table VIII. Note that the results in Table VIII are also used in Ref. [45]. 
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We use the parameters c and d in Table VIII to determine the normalized slope 
— (5 In//c? In m2) for each ensemble. Although the dependence of the normalized slopes 
on the light daughter-quark mass is resolvable, the effects are small for the ensembles we 
use in the analysis (with light daughter-quark masses ranging from 0.05ms to 0.4ms). We 
expect similarly small effects from the spectator-quark masses. We also expect that the 
lattice-spacing dependence of the normalized slopes is small, because it is a dimensionless 
ratio. We therefore correct each lattice form factor in each ensemble by a factor 


1 + 


/ 


d\n f / m2 
5 In m2 V^2 


(2.30) 


where m2 and m2 are the kinetic masses corresponding to the simulation and tuned k^, 
respectively. The resulting relative shift for each ensemble is shown in Table VII. Although 
the corrections to Kb itself are signihcant for some ensembles, the corresponding corrections 
to the form factors are much smaller (< 2.3%), as a consequence of the small normalized 
slopes. 


III. CHIRAL-CONTINUUM EXTRAPOLATION 

Here we extrapolate the form factors at four lattice spacings with several unphysical 
light-quark masses to the continuum limit and physical light-quark mass. We use heavy- 
meson rooted staggered chiral perturbation theory (HMrSxPT) [46, 47], in the hard-pion 
and SU(2) limits. We also incorporate heavy-quark discretization effects into the chiral- 
continuum extrapolation. 


Table VIII. Fitted normalized slopes to all momenta for and their combined fit with 

daughter-quark mass dependence. 



B - 

TT 

B - 

K 

Combined fit with rrii/ms 

dependence 


p value 

a In/ 
d In m2 

p value 

d\nf 
d In m2 

p value 

c 

d 

f± 

0.56 

0.145(11) 

0.46 

0.115(9) 

0.4 

-0.146(11) 

0.032(3) 

f\\ 

0.51 

0.160(16) 

0.74 

0.139(13) 

0.84 

-0.165(17) 

0.025(8) 

fr 

0.66 

0.146(16) 

0.79 

0.126(13) 

0.88 

-0.137(17) 

0.034(8) 
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A. SU(2) staggered chiral perturbation theory in the hard-pion limit 


The full-QCD next-to leading order (NLO) HMrSyPT expression for the semileptonic 
form factors can be written [46] 

yNLO ^ jjo) -t- c(xval + CaXsea + C^XE + c(Xe + ^iXa'^] , (3.1) 

where J =T, ||,T. Note that the expressions are in units of the mass-independent scale ri 
and the coefficients c/ have the dimension of The leading-order terms are 

f{0) _ _ 9B*Bn _ 

C = 4. (3.3) 

Jtt 

with gB*B-K the coupling constant and A^* = Mb* — Mb the B*-B mass splitting. 

The terms 5/j,iogs and SDi^gs are the one-loop nonanalytic contributions in the chiral ex¬ 
pansion, and depend upon the light pseudoscalar meson mass and energy [46]. Note that 
in the heavy-quark expansion Jt is proportional to f± up to 0{l/mb). We therefore use 
the same pole location and nonanalytic corrections for fx as f±. The terms analytic in Xi 
are introduced to cancel the scale dependence arising from the nonanalytic contribution in 
Eq. (3.1). The dimensionless variables Xi proportional to the quark mass, pion energy, 
and lattice spacing. We define 


2/im; 

Wal - g^2J2’ 

(3.4) 

_ /^(2m; + m^) 

!X^sea Q 9 pn ? 

SttV; 

(3.6) 

Xs = . . , and 

4vrA 

(3,6) 

AA 

“ STT^ff 

(3.7) 

Note that the valence mass mi is equal to the sea mass in our data. 

constant fi relates the pseudoscalar meson masses to the quark masses. 

The low-energy 

^|,ps = (”^i + ^ 2 ) 9 - + a^Ag, 

(3.8) 


and Ag is the mass splitting for staggered taste The average taste splitting in Eq. (3.7) 
is A = The quantities /r and A^ are obtained from the MILC Collaboration’s 

analysis of light pseudoscalar mesons and are shown in Table IX. 
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We constrain the parameter gB*B-K with a prior. The value of gB*B-K has been calculated 
with lattice QCD in the static limit [48, 49] or with a relativistic b quark [50] on gauge helds 
generated with domain-wall or Wilson sea quarks [51]. We set the prior, based on these 
lattice-QCD calculations, to be gB*B-K = 0.45 ± 0.08, where the error covers the differences 
among different determinations of the coupling. The LO and NLO coefficients, {cj, 0 < i < 
5}, are well determined by the data. Note that the formula given in Eq. (3.1) is slightly 
different from that in Ref. [5] where the NLO coefficients therein are our |c//cq| [i ^ 0). 
With the introduction of variables Xi defined in Eqs. (3.4)-(3.7), we should expect that 
Io^/cqI ^ 1 (* 7^ 0), or |c/| < |co|. In the actual hts, < 0.6 and < 1.0. 

Note that the coefficients cj are dimensionful, and they are evaluated here in ri units. We 
constrain them with loose priors: = 0 ± 1 and = 0 ± 2. 

Standard HMrSyPT uses the assumption that the external and loop pions are soft, he., 
[52, 53]. In our work, however, the external pion energies can be quite large, in 
some cases as much as 7 times the physical pion mass, and standard HMrSyPT may not 
converge well enough in this range. Indeed, the fit of the lattice form factor /y to Eq. (3.1) 
gives a poor conhdence level (p ~ 0), which is not improved by including higher-order con¬ 
tributions in the chiral expansion. Bijnens and Jemos [54] proposed an approach called 
hard-pion yPT, in which the internal energetic pions are integrated out and the depen¬ 
dence is absorbed into the low energy constants. ^ Since hard-pion yPT provides a more 
appropriate description of our data, we adopt it in this analysis. The explicit expressions 
for the hard-pion nonanalytic terms using SU(3) chiral perturbation theory as well 

as its SU(2) limit are given in the appendix of Ref. [45]. We take the SU(2) limit by inte¬ 
grating out the strange quark. The resulting expression has no explicit strange-quark mass 
dependence, which has been absorbed into the value of the low energy constants. The SU(2) 
hard-pion yPT provides a better description of our /y data than the SU(3) hard-pion yPT 
{p value 0.29 versus 0.09 from the NLO yPT ht with priors). We also hnd that the chiral 
expansion converges faster using SU(2) yPT when including higher-order chiral corrections 
in the ht to our data, which results in smaller yPT truncation errors than from using SU(3) 
yPT. Finally, Ref. [53] provides phenomenological arguments to prefer the application of 
SU(2) HMyPT over SU(3) to lattice-QCD data. We therefore use the SU(2) formula for our 

^ The factorization of hard-pion yPT breaks down starting at three loops [55] , but we only use the one-loop 
non-analytic terms. 
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Table IX. Fixed parameters that enter the chiral-continuum extrapolation fit function. The taste 
splittings rfa^Ag, ^ = P,A,T, V,I are for the pseudoscalar, axial-vector, tensor, vector and scalar 
tastes, respectively. The pseudoscalar taste splittings are zero by virtue of the remnant chiral 
symmetry of staggered fermions. The hairpin parameters were determined in a combined 

fit to light-pseudoscalar quantities at multiple lattice spacings. We take the result for these couplings 
at the a ~ 0.12 fm lattice and scale them to other lattice spacings by the ratio Ai.ins(o)/Arms(0.12fm) 
where Arms is the rooted mean square of the taste splittings. The continuum value of the low-energy 
constant ^ is evaluated at the same scale as the 0.09 fm lattice in our mass-independent scheme. 


a(fm) 0.12 0.09 0.06 0.045 

0 

rfo^Ap 0 0 0 0 

0 

rlo?AA 0.22705 0.07469 0.02635 0.01041 

0 

rfa^Ap 0.36616 0.12378 0.04298 0.01698 

0 

rfa^Ay 0.48026 0.15932 0.05744 0.22692 

0 

r^a^A/ 0.60082 0.22065 0.07039 0.02781 

0 

rla?6'y 0.0 0.0 0.0 0.0 

0 

rla^S'^ -0.28 -0.09 -0.03 -0.01 

0 

rip 6.83190 6.63856 6.48665 6.41743 

6.015349 


central value fit, but also check the consistency with the SU(3) hts in Sec. IV. 

Based on the above discussion, we use the following conditions for f±, /y and Jt in 
Eq. (3.1): 

j- £ _ jr £ hard, SU(2) 

v/j,logs — O/yiogs ’ 

0 , 

4 = 

2{2fimi) a‘^Aj/3 
^vai = g^2j2 “ SttV^ ’ 

where Eq. (3.10) is a consequence of the hard-pion limit, Eq. (3.11) and the factor 2 in 
the first term of Eq. (3.12) follow from the fact that we take mi = m'l and m'f^ has been 
integrated out, Eq. (3.12) preserves the chiral scale independence of the SU(2) hard-pion 
NLO expression, and a^A/ is the taste splitting of the taste-singlet pseudoscalar meson mass. 
The hts of the lattice form factors using NLO SU(2) hard-pion HMrSxPT have acceptable 


(3.9) 

(3.10) 

(3.11) 

(3.12) 
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confidence levels. We find, however, that there is a sizable shift in the £t result when 
including higher-order terms in the xPT expansion. We therefore need to study the effects 
of higher-order contributions in the chiral expansion. 

B. Next-to-next-to-leading order (NNLO) corrections 

We supplement the NLO SU(2) hard-pion yPT expression with the following NNLO 
analytic terms 

<^/Satytic = CeXvalXE + C'^jXa^XE + C^x| + 

+CioXvalx| + CiiXa2Xval + c(2Xa^XE + ^ 13 X 02 + ^mX^j (3.13) 

such that the complete NNLO xPT expression is, 

= /r° + (3Ai) 

Note that here uses the hard-pion and SU(2) xPT) manifested in Eqs. (3.9)-(3.12). 
All light-quark discretization errors that arise from taste violations are included here; generic 
errors from light-quark and gluon action, which are 0{asa?^^), are discussed in Sec. IIIC. 

Again, the expectation from chiral perturbation theory is that the coefEcients of these 
analytic terms should satisfy |c//cq| ~ 0{1) when written in terms of the dimensionless 
variables x given in Eqs. (3.4)-(3.7). We set the priors for the NNLO coefEcients for /±,||, 
{Cj ’,6 < i < 14}, to be 0 ± 0.6, since the role of these terms is simply to absorb the 
efiects of higher-order contributions in the chiral expansion. This width 0.6 corresponds 
to the size of |cq’V^^^|. For the same reason, we set the priors for the coefficients for /t, 
{cfr^^^,6 < i < 14}, to be 0 ± 1.0. Doubling the prior widths leads to negligible shifts on 
the central values of the form factors and less than 20% increases in the fit errors. 

C. Heavy-quark discretization effects 

The chiral-continuum extrapolation implemented in Eq. (3.14) accounts for the discretiza¬ 
tion effects from the gluons and the light staggered quarks. Discretization effects from the 
heavy b quark need a separate treatment. Heavy-quark discretization errors arise from the 
short-distance mismatch of higher-dimension Lagrangian and current operators [41, 42]. By 
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power counting, such mismatches are of 0{a?K^) or 0{asaA) where A is a QCD scale ap¬ 
propriate for the heavy-quark expansion. We follow the same method for incorporating 
the heavy-quark discretization effects described in Ref. [37] and include the following error 
function in Eq. (3.1), 

'^/j^ = {^eIe + + z^fy) {aAf + {zifs + ^ 3 / 3 ) (a^aA) + (3-15) 

where the mismatch functions fE,x,Y,B ,3 are given in the Appendix of Ref. [37]. The error 
functions Jb^e arise from mismatches of operators in the Lagrangian, while functions fx,Y,3 
arise from those of the vector current. The last term in Eq. (3.15) accounts for higher order 
heavy-quark and generic light-quark and gluon errors not included in Eq. (3.14), which is of 
the order ^^(aA)^. The £t parameters are constrained with priors: Oil for zy,ZB,Zo and 
0 ± \/2 for zx , zs'f the latter two are wider because the functions fx and /s both appear 
twice [42]. 

To summarize, after incorporating the heavy-quark discretization effects, the complete 
NNLO SU(2) hard-pion HMrSxPT expression is 

^NNLO+HQ ^ ^NNLO ^ + 6 ff^) , (3.16) 

where jg dehned in Eq. (3.14). With this treatment, the uncertainty due to truncating 

the chiral expansion at NNLO (cf. Sec. IV below), NNLO light-quark and gluon discretiza¬ 
tion effects, and LO heavy-quark discretization effects are incorporated in the £t error of 
the chiral-continuum extrapolation. The hts for f±, /y, and fy to Eq. (3.16) are shown in 
Fig. 11. 

To examine the size of discretization effects, we plot the form factors f± and /y with 
light-quark mass at each lattice spacing versus in Fig. 12. As we can see 

from the plots, the observed lattice-spacing dependence is very mild, with the data points 
at the largest lattice spacing (a ~ 0.12 fm) only about two statistical sigma away from the 
continuum limit. 

IV. SYSTEMATIC ERROR BUDGET 

The error output from the central-value ht described in Sec. Ill C already includes the 
systematic errors due to the light- and heavy-quark discretization effects and the uncertainty 


30 


X^/dof= 34.9/36, p = 0.52 



a~0.12 fm O.lOrrij 
a~0.12 fm 0.14mg 
a~0.12 fm 0.20mg 
a=0.09fm O.OSm^ 
a=0.09fm O.lOm^ 
a=0.09fm O.lSm^ 
a=0.09fm 0.20mj 
a=0.06fm O.lOm^ 
a=0.06fm Q.Mm^ 
a=0.06fm 0.20mj 
a=0.06fm 0.40mj 
a~0.045fm 0.20mg 
cont. phys. limit 


X^/dof= 52.8/48, p = 0.37 


X^/dof= 34.6/36, p = 0.53 




Figure 11. Chiral-continuum extrapolation of lattice form factors f± (upper left), /y (lower left) 
and fx (lower right) as functions of E.,^, where all quantities are in ri units. The colors denote the 
lattice spacings: 0.12 fm (gold), 0.09 fm (green), 0.06 fm (blue) and 0.045 fm (violet). The symbols 
denote the light-quark masses 0.05 (diamond), 0.1 (circle), 0.15 (square), 0.2 (downward¬ 

pointing triangle), and 0.4 (upward-pointing triangle). The colored lines correspond to the fit results 
evaluated at the parameters of the ensembles. The physical-mass continuum-limit curve is shown 
as a black curve with cyan error band. 


on gB*BTT- We now discuss other sources of systematic uncertainty. We tabulate systematic 
error budgets for /_|_ and /o at a representative kinematic point = 20 GeV^ within the 
range of lattice data in Table X. We also present the error budget for the full simulated 
lattice momentum range in Fig. 17. 
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Figure 12. Discretization effects in the form factor f± (left) and /y (right) at a few kinematic 
points. The plots show the form factors on the m[ = 0.2m)j ensembles at each lattice spacing vs. 
for various pion momenta (a slight extrapolation/interpolation is applied to adjust the raw data to 
the same Ej^ri). The range E.,^ri G [0.28,1.2] is used in the extrapolation to the full kinematic 
range. 

A. Chiral-continuum extrapolation 

As discussed above, our central fit uses NNLO^ SU(2) hard-pion HMrSxPT including 
contributions from heavy-quark discretization effects and the uncertainty in Here 

we consider variations of the ht function and the data included to estimate truncation and 
other systematic effects. 

First, we study the effects of truncating the chiral expansion by adding next-to-NNLO 
(NNNLO) analytic terms ^/^^aiyUc ^^s with coefficients constrained with the same 

priors as the NNLO coefficients. The variations in /+ due to changing the order of the 
yPT analytic terms are shown in Fig. 13. The hts of different orders are consistent in the 

region where most of the simulation data are located. Although the central values and 
errors differ noticeably between the NLO and NNLO hts, the central values and errors of 
the NNNLO ht are very close to the NNLO ht, indicating that the chiral extrapolation has 
stabilized by NNLO. As discussed earlier, the NLO coefficients are well determined by the 
data and we use well-motivated priors based on expectations from xPT for the NNLO and 
higher order terms. The fact that the error saturates with NNLO shows that the preferred 

^ NLO + NNLO analytic terms. 
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Figure 13. The function (1 — q'^/M'^t)f+ obtained from the chiral-continuum fits with NLO log¬ 
arithms and analytic terms through NLO (hatched black band), NNLO (solid cyan band) and 
NNNLO (dashed red lines). Errors shown are fit statistical errors only. The vertical line indicates 
the lowest of the lattice data. 

fit already incorporates the uncertainty from truncating the chiral expansion, and that we 
do not need to add an additional systematic error. The NNLO £t error as a function of 
is shown in Fig. 14. 

The standard soft-pion HMrSyPT fits of f± have reasonable conhdence levels, but those 
of /|| do not. Here we estimate the effect of using the hard-pion formalism by using standard 
HMrSyPT for fj_ but still employing hard-pion xPT for /y. The resulting difference from 
the preferred fit is small, less than 1% for /+. The same conclusion also holds for the form 
factor /o. 

We use SU(2) yPT, instead of SU(3) yPT, for our central £t. To estimate the effect of 
this choice, we restore the strange-quark dependence of the logarithm and analytic terms in 
Eq. (3.1). A practical issue arises with NNLO SU(3) yPT, where the terms proportional to 
the sea-quark mass, Xsea, are not well constrained by our data because the strange sea-quark 
mass is so similar on all of our ensembles. To obtain some sensitivity to Xsea, we include 
data on an additional a ~ 0.12 fm ensemble with an unphysically small strange-quark mass, 
am'h gea = am'i = 0.005. With the inclusion of this ensemble, we End the £t parameters for 
the terms involving ygea are better constrained. The differences between the NNLO SU(3) 
hts and the preferred hts are shown in Fig. 14. For /+, the difference is within the statistical 
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Figure 14. Variations in the chiral-continuum extrapolation from different fit Ansatze. The shaded 
area shows the fit error from the preferred NNLO SU(2) fit. The other curves show the systematic 
deviations from the NNLO SU(2) fit under the variations discussed in the text. 

error. For /o, the difference lies outside the statistical error for some of the simulated 
range, but the NNLO SU(3) fit quality is poor, with a p-value of 0.01. Because SU(3) yPT 
does not provide a good description of our data for /o, we do not take the difference between 
NNLO SU(2) and SU(3) fits as a systematic error. 

To check how our results are affected by data with high momenta, we also perform a 
fit excluding data with p = (27r/L)(l, 1,1). As shown in Fig. 14, the form factors /+ and 
/o from the low-momentum fit agree very well with those from the preferred full-data fit 
for the region > 20 GeV^. The systematic difference increases for small q^, where the 
highest-momentum data provide important information. 

Figure 14 summarizes the effects of all these variations. Comparing the deviations be¬ 
tween the central values of the alternate and preferred fits to the statistical error of the 
preferred fit, we find that the deviations are almost always smaller than the statistical error 
of our preferred fit. This confirms that fit errors of our preferred fits adequately account 
for the systematic effects associated with these variations. We therefore do not quote any 
additional systematic error due to these sources. 

We include heavy-quark discretization effects in our chiral-continuum extrapolation. As 
a consistency check, we compare our result with a power counting estimate obtained by 
evaluating in Eq. (3.15) at the a ~ 0.045 fm lattice spacing, setting the coefficients 
Zi = 1 and taking A = 500 MeV for the heavy-quark scale. We find 6fj^ ~ 1.5%. Figure 
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Figure 15. Heavy-quark discretization effects in the chiral-continuum fit. The solid blue, solid 
black, and dashed red curves show the NNLO fit error, the preferred fit error, and the NNLO fit 
error added to a 1.5% power-counting estimate of heavy-quark discretization errors, respectively. 

15 shows that the NNLO ht error (without the heavy-quark discretization effects) added 
to the 1.5% power-counting estimate in quadrature yields a similar error to that of the full 
fit. Thus, again, it is not necessary to add an additional error to that of the preferred 
chiral-continuum ht. 

B. Light- and bottom-quark mass uncertainties 

The effect of mistuning the 5-quark mass in our simulation has been largely reduced via 
the corrections described in Sec. IIF. Errors still arise, however, from the uncertainty in the 
tuned value Kb itself and from the procedure for shifting the form factors. From Eq. (2.30) 
we estimate the relative error by 

^ ^ <^(1^2) ^ f d\nf \ l/m2-l/m2 

f \d\nfh 2 / l/m 2 \d\nfh 2 / l/m 2 ’ 

where 6 {l/m 2 ) is related to the uncertainty due to the error in Kb while 5(^p()^) is the 
uncertainty on the normalized slope. The values of the physical Kb with errors are given 
in Table VII, and we can hnd the statistical uncertainty of the normalized slope using 
Table VIII. Using Eq. (2.30), we hnd that the value of 6 f/f on all ensembles is at most 
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Figure 16. Subdominant systematic errors over the range of simulated lattice momenta. Error 
estimates are described in the text. 


0.6%. We take the average value for 6 f/f on all ensembles, which is 0.4%, to be the error 
due to tuning kj,, and assign the same error to /+ and %. 

To obtain the physical form factors, we evaluate the result of the chiral-continuum £t 
at the physical light- and strange-quark masses determined from the MILC Collaboration’s 
analysis of light pseudoscalar mesons [19]. (Although we use SU(2) xPT, we include an 
analytic term proportional to Xsea to allow for a slight shift to the physical strange sea- 
quark mass.) The errors on the physical mi = {mu + md)/^ and mg are 3.5% and 3.0%, 
respectively. We vary the light- and strange-quark masses at which the chiral-continuum fit 
function is evaluated by plus and minus one standard deviation, and hnd that it produces 
differences below 0.4% in both form factors. 


C. Lattice scale ri 

We convert the lattice form factors and pion energies to physical units using the relative 
scale Ti/a determined from the static-quark potential (see Table III) and the absolute scale 
ri = 0.3117(22) fm [37]. The statistical uncertainties on ri/a are negligible. We propagate 
the uncertainty in ri by shifting it ilex and repeating the chiral-continuum £t. We hnd 
shifts of at most 0.7% in the range of simulated momenta. 
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D. Current renormalization 


With the mostly nonperturbative renormalization procedure that we use for the heavy- 
light currents, there are two sources of error. The hrst is due to the nonperturbatively 
calculated flavor diagonal factors Zy4^ and Zy4. Their values and errors are given in Table VI. 
We estimate the systematic error due to the uncertainties of Zyi^ and Zyi by varying their 
values by one sigma and looking for the maximum deviations in the form factors /_|_ and /q. 
The resulting deviations are small, ranging from 0.4% to 0.6%. 

The second source of error is due to the truncation of the perturbative expansion in the 
calculation of the pj. Because the pj are dehned from ratios of renormalization constants, 
their perturbative corrections are small by construction. Indeed, as seen in Table VI, for 
V^i they are less than 1% and for V^i they range between 2-3%. For the scale-independent 
vector current, we observe that the one-loop corrections to pyA are smaller than those for pyi^, 
and we use the same error estimate for both. In order to accommodate possible accidental 
cancellations, we take the error as 2pmLxag, where pmLxtts is an upper bound of the one-loop 
correction to in the range of heavy-quark mass arriQ < 3 that corresponds to the range 
of lattice spacings included in our analysis. The coupling is evaluated at the scale of the 
next-to-hnest lattice spacing in our calculation, a ~ 0.06 fm. This procedure yields an error 
estimate of 1%, which is larger than the one-loop correction to py4 over most of the mass 
range, and amounts to about 50% of the one-loop correction to pyi^ in the mass range that 
corresponds to the three hnest lattice spacings. This leads to an error of 1% for both fy 
and /o due to the perturbative renormalization factors. 


E. Finite volume effects 

We estimate the size of the hnite-volume effects by replacing the inhnite-volume chiral 
logarithms with discrete sums and repeating the chiral-continuum extrapolation. The change 
in our preferred ht after including hnite-volume corrections is very small, less than 0.01%, 
which we simply neglect. 
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F. Summary 


Figures 16 and 17 visually summarize the systematic error budget for the vector and 
scalar form factors /+, /o in the simulated lattice-QCD momentum range. By far the largest 
contribution to the total uncertainty is from the £t error, which includes the statistical 
uncertainty in addition to the chiral-continuum extrapolation and heavy-quark discretization 
errors. The total error on /_|_ is smallest, about 3%, in the region of ~ 20-24 GeV^. 

The subdominant errors, such as those from heavy-quark mass tuning, the current renor¬ 
malization etc., have mild dependence, as can be seen in Fig. 16. We therefore treat them 
as constant in g^ when propagating them. For each source, we take the maximum estimated 
error in the simulated g^ range; we then add these individual error estimates in quadrature 
to obtain an overall additional systematic error 6f. We find 6f^ = 1.4% and 6fg = 1.5%. 

In the next section, we will use our result for /_|_ to obtain \Vub\ via a combined £t with 
experimental data to the 2 ; expansion. Due to phase-space suppression, the experiments 
have poor access to the large-g^ region. On the other hand, the lattice-QCD form factor 
has a larger error than experiment at small g^ due to the sizable g^ extrapolation. As 
discussed below, the value of \Vub\ is mostly determined in the region g^ k, 20 GeV^, which 
is at the low end of the g^ range where the lattice-QCD form-factor error is still small. We 
therefore provide tabulated error budgets for the two form factors /+, /o from our calculation 
at the particular kinematic point g^ = 20 GeV^ in Table X. The error on /+(20 GeV^) is 
approximately 3.4%, which is about one third of the error on our previously-determined 
form factor in Ref. [5]. 

We compare our results for /_|_ and % with full errors, which are obtained by adding the 
fit errors from the yPT hts and <5/ in quadrature, with previous lattice-QCD calculations 
in Fig. 18. Our result for /+ agrees with previous results obtained at g^ > 17 GeV^ from 
Refs. [4, 5, 13], but is more precise. Our result for % is consistent with Ref. [13], but not 
with Ref. [4]. 

V. z EXPANSION AND DETERMINATION OF \Vub\ 

The chiral-continuum extrapolation described in the previous sections yields the form 
factors in the range 17 GeV^ < g^ < 26 GeV^. In this section, we extrapolate them to the 
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Table X. Error budgets of form factors /+ and /q at = 20GeV^. 


Uncertainty 

5/+ 

< 5/0 

Statistical+xPT+HQ+ 5 (B* Bn 

3.1 

3.8 

Scale ri 

0.5 

0.7 

Non-perturbative 

0.4 

0.6 

Non-perturbative Zyi 

0.4 

0.4 

Perturbative p 

1.0 

1.0 

Heavy-quark mass mistuning 

0.4 

0.4 

Light-quark mass tuning 

0.4 

0.2 

Total 

3.4 

4.1 



17 18 19 20 21 22 23 24 25 26 17 18 19 20 21 22 23 24 25 26 

(GeV") (GeV^) 


Figure 17. The distribution of the errors for /+ (left) and /q (right) as a function of The 
different bands in the plot show the contribution of the error source to the sum of squared errors 
(left y axis). The corresponding error can be read off from the right y axis. 

full kinematic range using the model-independent z expansion. The form factors resulting 
from the chiral-continuum extrapolation are functions specified by a set of parameters. One 
could, in principle, incorporate the expansion with the yPT expansion from the outset 
(see, e.g., Ref. [56]). With such an approach, however, the coefficients of the 2 ; expansion will 
have a nontrivial dependence on mi and a that must be derived from the underlying chiral 
effective theory. Because the dependence of the coefficients on a and mi is unknown, we 
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Figure 18. Comparison of /+ (left) and /q (right) from this work with previous lattice-QCD 
calculations by HPQCD [4], Fermilab/MILC[5] and RBC/UKQCD [13]. 


instead carry out the extrapolation in two steps, taking the chiral-continuum extrapolated 
results and feeding them into the z expansion. We introduce a functional method to perform 
the 2 ; expansion. We also apply the 2 ; expansion to the experimental data and, after verifying 
that the fits to experiment and to lattice QCD are consistent, we carry out a combined ht to 
obtain | 14 fe|- A byproduct of the last step is a precise determination for /+(q'^) constrained 
by lattice QCD at high and experiment at low q^. 


A. 2; expansions of heavy-light semileptonic form factors 


The z expansion involves mapping the variable to a new variable z by [57] 

\/t+ - - \/4 - to 


z{t,to) = 




(5.1) 


where t± = {Mb ± and to is chosen for convenience below. This change of variables 
maps the whole complex plane onto the unit disk in the z plane, where the upper (lower) 
path along the branch cut [f+, 00 ) is mapped to the lower (upper) half of the circle enclosing 
the unit disk in the complex 2 ; plane. Choosing to = {Mb -k M^){^yMB — fyilfy-)^ centers 
the full kinematic range for semileptonic B —)■ niu decay around the origin z = 0, and, 
moreover, restricts z to \z\ < 0.28. The small, bounded interval, together with a constraint 
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Table XL The BCL constants used to estimate T,(b,Nz). 


Boo Boi Bo 2 Bo 3 Boi -Bos -Boe 
fo 0.1032 0.0408 -0.0357 -0.0394 -0.0195 -0.0055 -0.0004 
/+ 0.0198 0.0042 -0.0109 -0.0059 -0.0002 0.0012 0.0011 


from unitarity ensures convergence of the expansion. As discussed below, we find in practice 
that the convergence is rapid. 

The form factors /+ and fo are analytic in ^ except for the branch cut [t+, cx)) and poles 
in [t_,t+]. We can write 


Pi{z)(j)i{z)fi = 


= > (5.2) 

n 

where Pfz), i = +, 0, are the Blaschke factors, which are introduced to remove the poles of 
fi in the region and (j)i{z) are the outer functions [57, 58]. We choose simple outer 

functions 0+ 0 = 1 and employ the following formulas to expand the form factors 

Nz-l 


f+iz) = 


1 - q^{z)/Ml, 


Eh 




N, 


w) = 


0 n 


(5.3) 

(6.4) 


n=0 


Equation (5.3) is known as the Bourrely-Caprini-Lellouch (BCL) expansion [59], which is 
constructed to reproduce the threshold behavior at and the asymptotic behavior as 

—)■ ±cx). Equation (5.4) is a simple series expansion of fo in z. 

The BCL coefficients in Eq. (5.3) and (5.4) obey the unitarity constraint [57, 59] 


N, 

S(6, W) = ^ Bmnbmbn < 1, 

m,n=0 


( 5 . 5 ) 


where the element Bmn satishes Bnm = Bmn = -Bo|m-n| and depends on the choice of to [59]. 
We tabulate the values of Bok for the form factors /+, fo in Table XL The inequality saturates 
when Nz —)■ cx). Although we do not incorporate this constraint into our hts, we check that 
our results satisfy it. 


B. Functional method for the z expansion 

In previous work, we have used synthetic data points generated from the yPT £t as 
inputs to the z fit [5], but here we take a new approach. We exploit the facts that the 
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xPT expansion is linear in the fit parameters and that it contains only a hnite number of 
independent functions (see Eq. (3.1)). We construct a covariance function iP(zi, Z 2 ), dehned 
as the covariance of any pair of points (zi,Z 2 ), using the set of functionals from the yPT 
expansion. Our new approach is to formulate the z expansion using the eigenfunctions of 
an integral operator dehned from K(zi,Z 2 )- 

Let us start with the NLO xPT expression Eq. (3.1), as an example. Because f± and /y 
are linear in their coefficients and cf, we can express them both in the compact form 

fj{mi,ms,a‘^,E^) = Cj ■ Xj, (5.6) 


where 


Cj^ 


Co Cl C2 C3 C4 C5 



(5.7) 


A',; = 




(1 ^/jjlogs) Xval Xsea Xe xl 



(5.8) 


and where J =_L, || and the variables are dehned in Eqs. (3.4)-(3.6). Any linear combination 
of f± and /|| can be written as 


f = [^ v] 


f± 

/|| 



^^11 


(5.9) 


with ^,7] functions of The uncertainty of the function / is encoded in the uncertainty 
in the coefficient vector Cj. In all these expressions, we are only interested in the terms 
with (or g^) dependence and, hence, 2 ; dependence. We can now dehne the covariance 
function K{z,z') in some valid domain [zi,Z 2 \. Explicitly, 


K{z,z') = Y{zf -Coy-Y iz'), 


(5.10) 


where 

and Cov is the covariance matrix of the involved coefficients c;( 


(5.11) 


CoVrnn 


(5.12) 
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The covariance function K{z,z') is a Mercer kernel [60], and Mercer’s theorem ensures 
that there exists a set of orthonormal functions "01 (^) dehned over the domain [^ 1 ,^ 2 ], such 
that 


z) = Y,i \ilJi{z)'ilJi{z'), (5.13) 

where Aj, xjji are the eigenvalues and eigenfunctions of the operator induced by the 
integral equation, 

rz2 

Lx'ipiz) = / K{z,z')'ilj{z')dz'. (5.14) 

J z\ 

The form factor f{z) can naturally be expanded in the basis of iJi{z)'- we only need to 
project the expansions in Eqs. (5.3) and (5.4) onto the same basis. The process of finding 
the expansion coefficients is equivalent to minimizing the following function (in analogy 
to the usual function, replacing the sum over discrete points with an integral over a 
continuous variable): 

xL = [ dz [ dz' [p^^{z) - gf{b,z)] K-\z,z') [f^^^{z') - gf{b,z')] 

2 


N, 




IJ Zl 


2 = 1 


A,; 


dz[f^^^{z) - gf{b,z)]i;i{z) 

Wz —1 

e{^{z)ipi{z)d. 

_n J Z^ 


n=0 


where 




is the form factor function from the yPT fit expanded in terms of 'ipi, and 

Vz-l 

atip,z) = 


(5.15) 


(5.16) 


(6.17) 


n=0 


are functions rewritten from the functions defined in Eqs. (5.3) and (5.4). For brevity we 
dehne 


(^) = 


_ ( 1 \n-Nz ^ 


(5.18) 

(5.19) 


1 - q^{z)/Ml^ L 

To summarize, we expand any form factor function obtained from the chiral- 

continuum extrapolation in the basis formed by the eigenfunctions of its covariance function 
K{z,z'). We then project the z expansion onto the same basis. Finally, we solve for the 
expansion parameters bn by minimizing the function Xiat defined in Eq. (5.15). 
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C. Details on z expansion of the form factors 


In addition to the fit errors from the chiral-continunm £t, we also need to propagate the 
subdominant errors, which have very mild dependence. We treat them as constant in 
and add them in quadrature, obtaining = 1.4% and = 1.5%. To include this effective 
subdominant error to the £t, we slightly modify the covariance function dehned in Eq. (5.10) 
by 

K'{z, z') = Y{zf ■ Cov' ■ Y{z'), (5.20) 

where the new covariance matrix includes the subdominant error, 

“I” SjrCmCn- (b’21) 

In the function array Y(z) dehned in Eq. (5.11), only a relatively small number of the 
elements are independent functions. For example, there are 42 terms in the NNLO SU(2) 
yPT £t functions for (including the HQ discretization contributions). Many of them, 
however, are set to zero in the continuum limit or become constant once the light-quark 
mass is hxed at its physical value. In the end, the chiral-continuum extrapolated /+ is 
described by only 6 independent functions. For %, the number of independent functions is 
7. Although we work in the functional basis in which the covariance function K{z,z') is 
diagonalized, singular modes can arise because K{z, z') is built upon Covj, which itself may 
have singular modes. Figure 19 shows the spectra of the operator Lk for form factor /_|_ q. 
The spectrum of /o contains two very small eigenvalues ;< 10“^^, and they are well separated 
from the other modes. When we discard these two modes, the ht quality of the functional 
2 ; ht improves from p = 0.03 to p = 0.46. For /_|_, we do not need to apply any cut on the 
eigenvalues. 

We hrst consider separate hts of /_|_ and % without any constraints on the coefhcients of 
the 2 ; expansion. With W = 3, or three free parameters bo,bi,b 2 , we obtain a low confidence 
level, p = 0.05, for the ht to %. The analogous three-parameter ht for /+ results in an 
acceptable conhdence level, p = 0.3. With W = 4 we hnd good conhdence levels for both 
form factors as well as sizable changes in the central values and errors (compared to the 
Nz = 3 case). The results of unconstrained hts of /+ and % with several values of W for 
/o and /+ are given for comparison in Table XII. The kinematic constraint /+(g^ = 0) = 
/o(<?^ = 0) is satished automatically, as is shown in Fig. 20 (left). 
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Figure 19. Eigenspectrum of the kernel function K{z,z') for /o and /+. 




Z Z 

Figure 20. The 2 : expansion of the form factors /+ and /q without any constraints (left) and with 
the kinematic constraint. The fits use Nz = 4. 


The ^-expansion coefficients should approximately satisfy the unitarity bound Eq. (5.5). 
Figure 21 (right) shows the bootstrap-sample distribution of S(6,4) from the fits for /q. The 
unitarity condition is marginally satisfied. In the case of /+, the sum S(5,4) ~ 0.03 is much 
smaller than the unitarity bound. Reference [61] pointed out that the smallness is expected 
based on heavy-quark physics. The value of S(6,4) should be of order [A/rribY, which is 
about 0.013 with a conservative choice of A = 1 GeV. The bootstrap-sample distribution 
of S(6,4) from the fits for /_|_ is shown against the heavy-quark estimate in Fig. 21 (left). 
They are consistent with each other. 

The result of the fits of /+,o with the kinematic constraint are shown in Fig. 20 (right). 
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150 

100 

50 

0 

Figure 21. Bootstrap sample distribution defined in Eq. (5.5) for the z expansion coefficients of /+ 
(left) and /o (right), compared with the heavy-quark scaling estimate 0.013 (red vertical line in the 
left plot) and the unitarity bound (red vertical line in the right plot). The 1000 samples of fits use 
only the kinematic constraint and Nz = 4. 



z 


Figure 22. Effects of truncating the 2 ; expansion for /+. The areas indicated by the hatched band, 
colored shaded band, and the lines are results with Icj errors for Nz = 3,4, 5 2 ; fits, respectively. 

With this constraint, we again examine how the £t varies with higher order Nz- We hnd 
that the ht central valnes do not change signihcantly when we change W from 4 to 5, in 
contrast to the case from 3 to 4, as is shown in Fig. 22 and Table XIII. 
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Table XII. Results from unconstrained 2 : fits of /o and 



f+ 


/o 


N, 

3 4 

5 

3 

4 

xVdof 

1.1 1.0 

1.6 

3.1 

0.54 

dof 

3 2 

1 

2 

1 

P 

0.27 0.33 

0.2 

0.05 

0.46 

Brnn^m^n 

0.052(16) 0.024(39) 

2(4) 

0.60(24) 

2 . 0 ( 11 ) 

/(o) 

0.12(6) 0.23(20) 

0.40(34) 

-0.14(9) 

0.20(17) 

bo 

0.409(16) 0.409(12) 

0.407(15) 

0.493(22) 

0.510(23) 

hi 

-0.72(13) -0.63(20) 

-0.60(22) 

- 2 . 1 ( 2 ) 

-1.7(2) 

b2 

-1.0(2) -0.3(1.3) 

0.3(1.7) 

-0.8(5) 

1.2(9) 

bs 

1 ( 2 ) 

4(5) 


3(1) 

64 


7(7) 




We perform several additional checks to confirm the stability of our results against various 
£t choices. In our preferred fit, we set the integral range in Eq. (5.15) to be z = [—0.25, 0.01] 
(or equivalently = [19.8,26.0]). The results, however, do not change noticeably if we 
extend the integral range to z = [—0.249, 0.069] which covers the full range of simulated 
lattice momenta. This is because the statistical fluctuations and correlations of the form- 
factor functions are largely decided by the region —0.1 0, where the yPT fit results 

are the most precise. 

We also try removing the smallest eigenvalue from the covariance function K{z, z') for 
we hnd that the resulting central values are essentially unaffected. Finally, we also try the 
fit using, instead of the BCL formula, the Boyd-Grinstein-Lebed (BGL) formula, which uses 
more complicated outer functions [58]. We hnd that the resulting form factors are within 
one standard deviation of the BGL result. 

To summarize, we obtain our preferred result from a simultaneous fit to /+ and /o with 
W = 4 and with imposing the kinematic constraint /+(g^ = 0) = /o(g^ = 0). The results for 
the two form factors /_|_ and /o are plotted in Fig. 23. In this plot, the form factors obtained 
from the xPT ht are overlaid on the results of the z £t. The z £t faithfully reproduces the 
yPT hts in the region where xPT is reliable (indicated by the ranges of the hatched bands). 
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Table XIII. Results of simultaneous fits of /+ and /o with the kinematic constraint. The = 4 
fit is our preferred result. 


Fit 

Nz = 3 

X, = 4 

Nz = 5 

xVdof 

2.5 

0.64 

0.73 

dof 

6 

4 

2 

P 

0.02 

0.63 

0.48 

mn^ra n 

0 .11(2) 

0.016(5) 

1.0(2.3) 


0.33(8) 

2 .8(1.7) 

8(19) 

/(O) 

0.00(4) 

0.20(14) 

0.36(27) 


0.395(15) 

0.407(15) 

0.408(15) 

bt 

-0.93(11) 

-0.65(16) 

-0.60(21) 

bt 

-1.6(1) 

-0.5(9) 

-0.2(1.4) 



0.4(1.3) 

3(4) 

bt 



5(5) 

b^o 

0.515(19) 

0.507(22) 

0.511(24) 

b\ 

-1.84(10) 

-1.77(18) 

-1.69(22) 

hi 

-0.14(25) 

1.3(8) 

2 (1) 

bl 


4(1) 

7(5) 

bl 



3(9) 


The 2 ; coefficients with errors from our preferred fit and their correlation matrix are provided 
in Table XIV. This information is sufficient to reproduce the lattice form-factor results over 
the full kinematic range. 

Figure 24 shows a comparison of our results with other theoretical calculations of the 
form factors [13, 62]. While our results are consistent with the previous results, ours are 
significantly more precise in the region of z < 0.1. 

Finally, it is interesting to compare the lattice form factors with theoretical expectations 
from heavy-quark symmetry. In the soft-pion limit, the vector and scalar form factors /+ 
and /o are related as [63] 


hm 

/+(g2 
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Jb 

fs* 


gB*BTT 


(5.22) 









Table XIV. Central values, errors, and correlation matrix of the coefficients of /+ and /o from the 
Nz = 4 lattice-only z-fit with the kinematic constraint. 



K 

bt 

^3^ 

OO 

ft? 

50 


0.407(15) 

-0.65(16) 

-0.46(88) 0.4(1.3) 0.507(22) 

-1.77(18) 1.27(81) 4.2(1.4) 

6+ 1 

0.451 

0.161 

0.102 

0.331 

0.346 

0.292 

0.216 

K 

1 

0.757 

0.665 

0.430 

0.817 

0.854 

0.699 

K 


1 

0.988 

0.482 

0.847 

0.951 

0.795 

4 



1 

0.484 

0.833 

0.913 

0.714 

oo 




1 

0.447 

0.359 

0.189 

6? 





1 

0.827 

0.500 

too 






1 

0.838 








1 



z 


Figure 23. z-fit results for the form factors /o and /+ as functions of z. 

up to corrections of 0(1/ml). This expression updates the leading-order result of Ref. [64] 
to include the l/m;, correction, which turns out to be simply the additional multiplicative 
factor (/b*//s)~^ in fhe soft-pion limit. In Fig. 25 we plot the ratio of (/o//+)/(l — 
obtained using the coefficients of our preferred ^-expansion in Table XIV. We also show the 
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(1-q2/MB*2)f,(z) 


1.4 

1.2 

1 

_ 0.8 
0.6 

0.4 

0.2 

0 

- 0.2 - 0.1 0 0.1 0.2 
z 

Figure 24. Comparison of /+ (left) and /q (right) from the z-expansion fit results of this work with 
recent theoretical calculations using light-cone sum rules (LCSR) [62] and lattice QCD [13]. 

theoretical expectation from Eq. (5.22), taking the HPQCD Collaboration’s recent three- 
flavor lattice-QCD result for the decay-constant ratio /b*//b = 0.941(26) [65], and using the 
same value of gB*B-K = 0.45(8) as in our chiral-continuum extrapolation. The large width of 
the expected band is due to the generous range taken for Qb^b-k- Higher-order corrections 
in the heavy-quark expansion are expected to be small. Taking a conservative value for 
A = 500 MeV and rub = 4.2 GeV, one would estimate (A/rub)^ corrections to be about 1%. 
The difference of /b*//s from one also provides a measure of A/mb ~ 6%, which would 
indicate that (A/mb)^ corrections may even be below the percent level. The lattice form 
factors agree with the theoretical expectation for > 27 GeV^. 


T 


RBC/UKQCD2015 C 
This work 


- 0.2 - 0.1 


0 

z 


0.1 


0.2 



D. Determination of \Vub\ 

We now combine our lattice form factors with experimental data for B —)■ to obtain 

\Vub\- The Standard-Model partial branching fraction is TBdT/dq^, where dV/dq^ is defined 
in Eq. (1.1). The contribution from /o is negligible due to the small lepton mass. Given 
/+(g^), the branching fraction in the Ah q^ bin is 

AHf = cl\v^b? r" (5.23) 

J'il 
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Figure 25. The lattice form-factor ratio (/o//+)/(l — q^/M‘^t) versus (cyan shaded band) 
compared with the prediction in the soft-pion limit from heavy-quark symmetry and yPT [63] 
(hatched band). The width of the hatched band reflects only the uncertainty from gB*B-K = 0.45(8) 
and not other theoretical errors. 

where C% = (rBG|’)/(247r^) is a constant. For the combined lattice plus experiment z fit, 
we define a for fhe experimental measurements as 

xLp = - Ae«‘)Cov“yAJ5-'' - ABf), (5.24) 

where is the experimentally-measured branching fraction in the ith bin {i is a 

shorthand notation for each bin in each experiment included in the fit) and Cov®^^ is the 
experimental covariance matrix, including the statistical and all systematic errors. 

We use the experimental results compiled by the Heavy Flavor Averaging Group (HFAG) 
[6]; BaBar untagged 6-bin analysis (2011) [7], Belle untagged 13-bin analysis (2011) [9], 
BaBar untagged 12-bin analysis (2012) [8] and Belle tagged analysis with 13 bins for the 
and 7 bins for the B~ mode (2013) [10]. For convenience in the £t, we assume isospin 
symmetry to convert the Belle tagged B~ data to the B^ mode via 

AB(B° ^ vr+rn)Beiie b- = 2—AH(H~ ^ 7r°rz/), (5.25) 

tb- 

where r^o = 1.519(7) ps and Tb- = 1.641(8) ps are from the PDG [66]. 

We omit systematic correlations between the BaBar and Belle analyses, because they 
do not share any major systematic errors. The BaBar 6-bin and 12-bin data have very 
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small overlaps in the selection of samples, so the statistical errors can be considered ap¬ 
proximately uncorrelated. There is some systematic correlation between the two analyses, 
which is, however, supposed to be insignihcant [67]. The Belle untagged and tagged data 
are also largely uncorrelated because the dominant source of systematic errors in these two 
measurements are very different. In summary, we take the four experimental analyses as 
independent measurements. 

On the other hand, there are systematic correlations between the two isospin modes of 
the Belle tagged data, which we estimate as follows. Let and be the branching 
fractions in the ith and ath bin of the charged and neutral decay modes, respectively. 
Let a~ , (T° be the systematic uncertainties of the two modes from source x and rbe 
the correlation between them. Then we estimate the off-block-diagonal elements of the 
systematic error covariance matrix by 

Si.= E r-" , (5.26) 

xGall sys. 

where the sum is over all sources of systematic errors. That said, only a few of the systematic 
errors contribute noticeably to the sum and the biggest source of error, the tag calibration, 
dominates. From the correlation matrices, we construct the total covariance matrices of each 
isospin decay mode by adding the statistical matrices and the systematic matrices. We then 
take the direct sum of the covariance matrices of the B~ and modes block-diagonally 
and add the off-block-diagonal elements Sia so that we can £t them simultaneously. 

We hrst £t the 2 ; expansion to the experimental data only and without any constraints 
on the coefficients. We use the BCL formula with three parameters, W = 3, where the 
normalization is iWbl^o • The result is shown in Table XV. To check the consistency in the 
shape among the experimental data sets, we also £t each experimental data set separately. 
The individual hts all have acceptable conhdence levels and p values, but the combination of 
all four data sets gives a rather poor £t that is not improved by going to higher order in 
e.g., Nz = 4. The poor ht stems from the BaBarll measurement, which is only marginally 
consistent with the other three. Figure 26 compares the shapes (slopes 61/&0 and curvatures 
62 /& 0 ) of fhe separate and combined experimental hts with the lattice-only ht. The lattice 
form factor shape is consistent with all of the experimental results. 

To perform a combined ht to the lattice and experimental data, we dehne the total chi 
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Table XV. The results of fits to experimental data only. 


Fit 

xVdof 

dof 

P 61/60 62/60 

6 o|K 6| X 10-3 

All exp. 

1.5 

48 

0.02 -0.93(22) -1.54(65) 

1.53(4) 

BaBarll [7] 

2 

3 

0.12 -0.89(47) 0.5(1.5) 

1.36(7) 

BaBarl2 [8] 

1.2 

9 

0.31 -0.48(59) -3.2(1.7) 

1.54(9) 

Belle 11 [9] 

1.1 

10 

0.36 -1.21(33) -1.18(95) 

1.63(7) 

Bellel3 [10] 

1.2 

17 

0.23 -1.89(50) 1.4(1.6) 

1.56(8) 


• 

BaBarll 

• 

Bellell 

• 

BaBarl2 

• 

BellelB 

• 

All expt. 

0 

Lattice only 



Figure 26. The contour plot of the slope 61/60 and curvature 62/69 of the form factor /+. The 
open ellipses are the l-cr contour of the slope and curvature constructed from the 3-parameter z fit 
to individual experimental data. The gold filled ellipse is from the combined fit of all experimental 
data. The cyan filled ellipse is from the 4-parameter z fit to lattice form factors. 


squared function, 


2 _ 2,2 _ l 2 _ l 2 

X — Xlat + XsaBarll + XBellell + XBaBarl2 + XBellelS) 


(5.27) 


where the lattice and experimental chi squared functions are defined in Eqs. (5.15) and 
(5.24), respectively. The fit is performed to these five independent data sets with common 
shape parameters bm and overall normalization \Vub\ by minimizing Eq. (5.27). Table XVI 
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Table XVI. Results of the combined lattice+experiment fits with = 4;. 


Fit 

xVdof 

dof 

p value 

6 + bt bt bt 

|Wb|(xl03) 

Lattice+exp. (all) 

1.4 

54 

0.02 

0.419(13) -0.495(55) -0.43(14) 0.22(31) 

3.72(16) 

Lattice+BaBarll [7] 

1.1 

9 

0.38 

0.414(14) -0.488(73) -0.24(22) 1.33(44) 

3.36(21) 

Lattice+BaBaiT2 [8] 

1.1 

15 

0.34 

0.415(14) -0.551(72) -0.45(18) 0.27(41) 

3.97(22) 

Lattice+Bellell [9] 

0.9 

16 

0.55 

0.412(13) -0.574(65) -0.40(16) 0.38(36) 

4.03(21) 

Lattice+Bellel3 [10] 

1.0 

23 

0.42 

0.406(14) -0.623(73) -0.13(22) 0.92(45) 

3.81(25) 


summarizes the various fit results. Due to the tension between the experimental data sets, 
the p value of the fit to the lattice result and all experiments is only 0.02. Table XVII shows 
the contributions to the total from each data set of the combined fit. By far the largest 
contribution to x^/dof is from the BaBar 6-bin data set, similar to what we find for the 
experiment-only fits presented in Table XV. 

In the combined fit to lattice form factors and experimental data, the kinematic constraint 
between /+ and /o at = 0 is unimportant for the determination of \Vub\- This is because 
the experimental data constrain the shape at low g^. Removing the kinematic constraint 
from the combined fit and fitting only with the vector form factor /_|_ changes neither the 
coefficients of the z expansion nor the value of \Vub\- We also try varying the number of 
parameters bm in the expansion (W). The results are shown in Table XVIII. Compared 
to our preferred fit with W = 4, the fit using W = 3 gives a very low p value and a shift of 
about la in both the form factor and \Vub\-, while the fit result using = 5 nearly coincides 
with that of the W = 4 fit and the values of \Vub\ are almost identical. 

The experimental data are plotted in Fig. 27 (left) along with the z fits to the lattice 
data and to all experimental data. The lattice form factor and experimental measurements 
provide complementary information and, when combined, yield an accurate description of 
the form factor over the full-g^ range and hence a precise determination of \Vub\- The plot 
shows that the experimental data dominate the determination of the form-factor shape in 
the large-z (small-g^) region while the lattice-QCD form factor dominates the small-z (large- 
g^) region. In the intermediate region around g^ ~ 20 GeV^ {z ~ 0), the lattice-QCD and 
experimental uncertainties are similar in size. This region is decisive in determining \Vub\ 
and, hence, can be used to estimate the separate contributions from lattice and experimental 
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Table XVII. The contribution to the total from each data set of the combined fit. The column 
data” gives the number of independent functions (for lattice QCD) or the number of bins (for 
experiment). The total X^/(# data) agrees with the x^/dof in Table XVI, once the constraint and 
number of fit parameters have been taken into account. 


data set 

^ data 


x7# data 

Lattice 

11 

4.8 

0.44 

BaBarll [7] 

6 

20.9 

3.5 

BaBarl2 [8] 

12 

15.1 

1.3 

Bellell [9] 

13 

13.8 

1.1 

Bellel3 [10] 

20 

23.5 

1.2 

Total 

62 

78.2 

1.26 


Table XVIII. Combined lattice+experiments 2 ; fits with = 3,4 and 5. 


N, 

X^/dof 

dof 

p value 

6+ 6+ 6+ 6+ 

ht \Vub\ 

3 

2.5 

56 

0.0 

0.425(12) -0.424(31) -0.59(9) 

3.63(11) 

4 

1.4 

54 

0.02 

0.419(13) -0.495(55) -0.43(14) 0.22(31) 

3.72(16) 

5 

1.5 

52 

0.01 

0.418(13) -0.491(56) -0.31(30) 0.01(55) 

-0.6(1.9) 3.72(16) 


data to the \Vub\ uncertainty. At = 20 GeV^, the error on the lattice-QCD form factor /_|_ 
is about 3.4% (see Table X) and the error on f^\Vub\ from the experiment-only £t is 2.8% 
at the same momentum. Adding these two errors in quadrature gives a total uncertainty of 
4.4%, which is consistent with the error on \Vub\ obtained from the full £t, 4.3%. Another 
estimate of the individual error contribution to \Vub\ can be obtained from the uncertainty 
on the £t parameters from the separate lattice-QCD and experiment hts. From the £t to all 
experimental data in Table XV, the normalization is |Vu;,|6o = (1.53±0.04) x 10“^. Similarly, 
the lattice-only x £t gives the normalization % = 0.407 ± 0.015. Assuming no correlation, 
one would obtain |14fe| = (3.76 ± 0.17) x 10“^, which is close to what we obtain from the 
combined fit. 
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Figure 27. Left: comparison of vector form factor f+(z) from z expansion fits to: only the lattice- 
QCD data (cyan band) and only experimental data including all four measurements (gold band). 
Right: the similar plot for the partial branching fraction dB/dq^. The fits including lattice results 
use Nz = 4, while the experiment-only fit uses = 3. The experimental data points and the 
experiment-only z-fit result in the left plot have been converted from (Ai3/Ag^) to /+ using 
\Vub\ from the combined fit. The lattice-only fit result (cyan band) and the combined-fit result (red 
band) in the right plot is converted from the form factor with the same \Vub\- 

VI. RESULTS AND CONCLUSION 

Our final result for \ Vuh\, obtained from our preferred z fit combining our lattice-QCD cal¬ 
culation of the B —)■ 'K^v form factor with experimental measurements of the corresponding 
decay rate, is 

\Vuh\ = (3.72 ± 0.16) X 10"l (6.1) 

The error includes all experimental and lattice-QCD uncertainties. The contribution from 
lattice QCD to the total error is now comparable to that from experiment. The error reported 
here, following HFAG [6], does not apply the PDG prescription for discrepant data; that 
prescription [66] would scale the error by a factor of -y/x^/dof = 1.2. As can be seen from 
Table XVII and Fig. 26, the low £t quality is due to the tension between the BaBarll data 
set and the others. An inspection of all the experimental data in Fig. 27 shows that the 
point near 2 ; = —0.1 in the BaBarll data set is lower than the others and a bit more precise 
than one might have anticipated, but does not suggest that this or any of the data sets have 
any systematic problems. 
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We compare our determination of \Vub\ with other results in Fig. 28. In particular, our 
result is consistent with the recent determination from HFAG using our collaboration’s 2008 
form-factor determination [5] obtained from a small subset of the gauge-field ensembles used 
in this work. The difference in the central values is due to a small shift in the central values 
for the form factor /_|_ of this analysis compared to our previous analysis [5]. As shown 
in Fig. 18 (left), the form factor /_|_ from this analysis is consistent within errors with the 
previous analysis, but shifted slightly downward and with an error smaller by roughly a 
factor of three. The two analyses have very little statistical and systematic correlation. Our 
result is also compatible with Standard-Model expectations from CKM unitarity [70, 71]. 
Although our determination of \Vub\ is higher than that in Ref. [5], and thus closer to 
the determination from inclusive B —)■ semileptonic decays [6], the inclusive-exclusive 

disagreement is still greater than 2a. 

A byproduct of the combined lattice and experiment fit is a more precise determination 
of the vector and scalar form factors than from the lattice-QCD calculation alone. Both 
form factors /_|_ and /o are well determined from lattice QCD in the high region, and 
/+ is strongly constrained by experiment in the low region. This information is then 
transferred to /o via the kinematic constraint /o(0) = /+(0). The resulting form factors are 
shown in Fig. 29. The corresponding ^-expansion coefficients and their correlations are given 
in Table XIX. These represent the present best knowledge of the B —)■ form factors, and 

can be used in other phenomenological applications or to test other nonperturbative QCD 
calculations. 

Future improvements in the determination of the R —)■ tt semileptonic form factor /+ will 
further reduce the uncertainty on \Vub\- If the uncertainty of at ~ 20 GeV^ can 

be reduced further from 3.4% to 1.5%, we would expect a precision of 3% in |14fe|, using 
the current experimental input. With the anticipated improvement in the experimental 
rate measurement from Belle II, this error would be reduced further. The reduction of 
uncertainty in is expected with the newly-available MILC gauge ensembles that are 

being generated using the highly improved staggered quark (HISQ) action [72]. The new 
HISQ ensembles have statistics similar to the asqtad ensembles, but with much smaller 
light-quark discretization effects. Further, the HISQ ensembles simulated at the physical 
light-quark masses will remove the need for a chiral extrapolation, thereby eliminating a 
significant source of uncertainty in this work. These ensembles have already helped to 
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HPQCD 2006 + HFAG 2014, B nlv 
Detmold et al. 2015 + LHCb 2015, A^ ^ plv 
BLNP 2004 + HFAG 2014, B^Xlv 

U 

UTFit 2014, CKM unitarity 


Figure 28. Determinations of | Vub\- The squares are obtained from B —)• Triu decay using theoretical 
form factors from this analysis, our earlier work [5] (now superseded, but with updated experimental 
input from HFAG 2014 [6]), a three-flavor lattice calculation by RBC/UKQCD [13], light-cone sum 
rules (orange square) [62], and HPQCD [4] (using the > 16 GeV^ experimental data only). The 
blue upward-pointing triangle is obtained from A;, —>■ pia decay using lattice-QCD form factors 
from Ref. [68] and experimental data from LHCb [69]. The black diamond shows the inclusive 
determination using B — Xu^v decays [6] with the theoretical approach of Ref. [15]. Also shown 
is the expectation from CKM unitarity [70] (green filled circle). For the exclusive determinations 
from B —)• nil' decay (squares), all four experimental results [7-10] are used except in the LCSR 
z-fit where only the more recent BaBar [8] and Belle [10] data are used. 

determine the form factor [73] and the leptonic decay constants and fx [74], 

and hence the relevant CKM matrix elements [Kisl, \^cd\ and |V/s|, with high precision. All 
of these improvements will further rehne and reduce the uncertainties in \ Vub\, and may also 
help to resolve the inclusive/exclusive puzzle. 
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Figure 29. The form factors f-\-{z) (left) and foiz) (right) from the 2 ; fits to the lattice data (cyan 
band), to all experiments (gold band), and to the lattice data plus experiments (red band). The 
fits including lattice results use = 4 and the kinematic constraint, while the experiment-only 
z fit uses Nz = 3. The experiment-only z-fit result has been converted from (A.K/Ag^)^^^ to /+ 
using \ Vub\ from the combined fit. 
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Table XIX. Central values, errors, and correlation matrix of the coefficients of /+ and /q from the 
Nz = 4 z-Rt combining lattice and all four experiments. 



\Vub\x 10^ 

bt 

bt 

bt 

^3^ 



6° 



3.72(16) 

0.419(13) 

-0.495(54) 

-0.43(13) 0.22(31) 0.510(19) 

-1.700(82) 1.53(19) 4.52(83) 

\Vub\ 

1 

-0.870 

-0.400 

0.453 

0.428 

-0.175 

-0.201 

-0.119 

-0.009 

K 


1 

0.140 

-0.455 

-0.342 

0.224 

0.174 

0.047 

-0.033 

K 



1 

-0.789 

-0.874 

-0.068 

0.142 

0.025 

-0.007 

bt 




1 

0.879 

-0.051 

-0.253 

0.098 

0.234 

bt 





1 

0.076 

0.038 

0.018 

-0.200 

OO 






1 

-0.043 

-0.604 

-0.388 

ft? 







1 

-0.408 

-0.758 









1 

0.457 

b^ 
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